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THEORETICAL AND EXPERIMENTAL STUDY OF A NEW METHOD 
FOR PREDICTION OF PROFILE DRAG OF AIRFOIL SECTIONS 


By S. H. Goradia and D. E. Li 1 ley 
Lockheed-Georgia Company 


SUMMARY 


This report describes theoretical and experimental studies which were con- 
ducted for the purpose of developing a new generalized method for the prediction 
of profile drag of sin gle componen t airfoil sections. This method aims at solu- 
tion for the flow irTtHe^ake from the airfoil trailing edge to the large 
distance in the downstream direction; the profile drag of the given airfoil 
section can then easily be obtained from the momentum balance once the shape of 
velocity profile at a large distance from the airfoil trailing edge has been 
computed. Computer program subroutines have been developed for the computation 
of the profile drag and flow in the airfoil wake on UNI VAC 1108 computer. The 
required inputs to the computer program consist of free stream conditions and 
the characteristics of the boundary layers at the airfoil trailing edge or at 
the point of incipient separation in the neighborhood of airfoil trailing edge. 
The method described in this report is quite generalized and hence can be ex- 
tended to the solution of the profile drag for mul ti -component airfoil sections. 

Experimental measurements of velocity profiles and static pressure dis- 
tribution in the wake of sharp trailing edge test airfoil were obtained by use 
of a pitot static tube. Shear stress profiles in the airfoil wake were computed 
by the indirect method. A hot-wire anemometer was used for the measurements of 
velocity profiles in the wake of the blunt trailing edge test airfoil. The 
above measurements were obtained for a 15 percent thick modified NACA four digit 
airfoil section which was modified near the leading edge for the purpose of 
increasing the maximum lift coefficient. 

A physical flow model for the wake flow was developed from the experimental 
data. This flow model is divided into regions and layers depending upon the 
characteristics of velocity profiles. Parameters for velocity profile similari- 
ty functions were developed from experimental measurements in different regions. 
Integral equations were derived for the solution of wake flow and these equations 
are coupled, ordinary nonlinear and nonhomogeneous differential equations. 
Parameters, such as shear and static pressures at the edges of various layers 
appear in the above equations as coupling terms. Functional representation of 
these parameters in terms of dependent variables, arranged in dimensionaless 
groups, was accomplished with the help of the principle of local dynamic 
simi larity. 

The numerical solution of the above equations was accomplished by the 
modified Euler method with repeated iterations. Results of computations of the 
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computer program subroutines for the flow in the wake and profile drag were 
compared with experimental data for three airfoil sections. Reasonable agree- 
ment was obtained between measured and calculated flow quantities in the wake 
and profile drag corresponding to the occurrence of incipient separation on the 
upper surface of airfoil. 


I . INTRODUCTION 


With the recent energy crisis and shortage of fuel, it has become impera- 
tive to search for ways of reducing the drag of aircraft operating in the low 
and _hJgh subsoni c speed regi me, for their economical operation during cruise. 

The reduction in the drag of an aircraft, and hence saving in the fuel consump- 
tion, can be accomplished by proper design of the wing section. The effective 
design of wing sections, however, requires an accurate and reliable method for 
predicting prof i le drag. The trade-offs between aerodynamic and structural 
characteristics of airfoils must consider profile drag as a merit factor. A 
prediction method, which is based on theoretically sound principles consistent 
with real flow phenomena, can be used to evaluate candidate wing sections with- 
out resorting to the elaborate, expensive and time consuming wind tunnel tests. 
Moreover, the development of airfoils for special appl icat ions, such as minimum 
drag at specified lift conditions, requires a method with general applicability. 
The development of such a method was aimed during this study contract with NASA, 
Lang ley . 

In order to estimate the drag of a body or to design a body for minimum 
drag, it is necessary to understand the origin and mechanism of the aerodynamic 
drag produced by a body. D'Alembert's paradox states that in an inviscid fluid 
a body can experience no drag. This can be proved relatively easily by use of 
the momentum theorem. Why then does a body experience drag in a real fluid? In 
the case of two-d i mens iona 1 bodies, e.g. airfoil section at a positive angle of 
attack, the induced drag associated with lift is nonexistent, but profile drag 
or parasite drag of such bodies have a finite value depending upon free stream 
conditions and geometrical shape of two-dimensional bodies. The above mentioned 
profile drag of a two-dimensional body is composed of two parts, namely, the skin 
friction drag and the form or pressure drag. These parts may be of equal magni- 
tude or the one may completely overshadow the other, depending upon the shape of 
the body. The skin friction drag is the result of the shearing stresses in the 
fluid as it passes over the surface of the body. The form drag results from the 
unbalance in normal pressure forces around the body due to the presence of 
attached or separated boundary layer around the surface of a body. These state- 
ments are exemplified by the use of Figure 1-1; this figure shows the cases of 
flat plate at 0° and 90° angles of attack and an airfoil at an angle of attack 
corresponding to the desired lift condition. 

In Figure 1-1 (a) the drag is entirely the result of skin friction, whereas 
in Figure l“l(b) it is entirely form or pressure drag. If some means could be 
used to prevent the separation of the flow at the edges of the plate shown in 
Figure I - 1 (b) , the drag could be reduced to zero. The total or profile drag for 
the airfoil section of Figure I — 1 (c ) depends upon geometrical contour of the 
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(b) Plat Plate at 90* 
Angle of Attack 


Plat Plate at 0* 
Angle of Attack 



(c) Airfoil at an Angle 
of Attack 



Presence of Skin 
Priction Brag only 



Case (b) — Presence of Pressure 
Brag only 

Case (c) — Both Skin Priction 
and Pressure Brag 
Constitute Total Brag 




FIGURE 1-1 - EXAMPLE SHOWING PRESENCE OP VARYING PROPORTION OP SKIN FRICTION 
AND PRESSURE BRAG CONSTITUENTS OF TOTAL PROFILE DRAG 


body and free stream conditions such as Reynolds number, angle of attack, and 
Mach number. Combination pf the above factors affect physical parameters, such 
as location and value "of pressure peak on the upper surface of airfoil, exis- 
tence of laminar boundary layer separation or transition, and conditions and 
growth of turbulent or laminar boundary layer on the upper and lower surfaces 
of the airfoil near its trailing edge. These physical parameters then determine 
the total profile drag of a given airfoil section and also the ratio of its con- 
stituents, namely, skin friction and pressure drag. The relation between 
pressure and friction drag of an axial ly-symmetric afterbody behind a long cir- 
cular cylinder is of interest both from the point of view of practical applica- 
tion in design and the theoretical studies. Available experimental data for 
the drag coefficient, Cq, for blunt and sharp afterbodies, which are shown in 
Figures 1-2 (a) and l-2(b), reveal an interesting phenomenon. Figure 1-2 (a) 
shows the plots of measured drag coefficients, Cp, versus Lr/Rq for blunt and 
sharp afterbodies; in this Figure Lr denotes the length of the afterbody and 
R 0 denotes the base radius of the afterbody. The experimental data shown in 
this figure reveal that blunt or thicker afterbody has a lower drag coefficient 
than the sharp afterbody. This result is quite surprising because it is known 
from available experimental data on airfoil sections (for example NACA four 
digit series) that thicker sections have higher values of profile drag than 
lower thickness ratio airfoil sections. The above-mentioned discrepancy of 
experimental results of Figure l-2(a) can be explained with the aid of detailed 
analysis of experimental results shown in Figure I —2 (b) . This figure shows the 
breakdown of the Jocal total drag into local integrated values of skin friction 
drag and pressure drag. Values of integrals 1 ^ and l 2 , which are shown plotted 
as a function of X/Lr in Figure I —2 (b) for both thicker and sharp afterbodies , 
represent local integrated values of skin friction and pressure drag, respec- 
tively. The algebraic sum of Ij and l 2 at particular X location on the after- 
body is equal to the total drag of afterbody up to that X location. The local 
integrated value of the pressure drag l 2 for blunt afterbody increases up to 
X/Lr^0.75 and then starts decreasing up to the tail end because of better 
pressure recovery; this better pressure recovery is due to the absence of flow 
separation. For the sharp afterbody, however, the value of 1 2 increases con- 
tinuously up to the tail end because of flow separation at X/Lr^0.8. Thus, 
the sum of lx and l 2 at the tail end, which represents the total drag co- 
efficient, is higher for the sharp afterbody than for the blunt afterbody. 

The total drag of an arbitrary airfoil section, at least in theory, can be 
computed from the knowledge of pressure distribution and airfoil geometry. From 
the knowledge of either experimental or calculated pressure distribution on the 
surface of the airfoil and the airfoil geometry, the pressure drag can be cal- 
culated from the Cp versus Z/C relationship of that airfoil, where Z/C is the 
nond imens ional ordinate of the given airfoil section, the determination of the 
algebraic sum of the familiar thrust and drag loops gives the resultant pressure 
drag. The boundary layer characteristics and hence the skin friction drag can 
be calculated for the given pressure distribution by the use of pertinent 
boundary layer theories. In actual practice, however, it is extremely difficult 
to obtain accurate drag level by the above mentioned which is defined as the 
Direct Method. Slight errors in pressure distribution, while not significantly 
affecting the boundary layer characteristics, lift and pitching moment, can 
create errors in the thrust and drag loops that are seriously magnified when 
calculating pressure drag by taking the difference of these two loops. Hence, 




= Local Integrated Value for Skin Friction Drag 
x/ l r L 

= / 2 JL JL m d 1*1 

° r q Ro q l r 


1- = Local Integrated Value for Pressure Drag 
x/L r 

= f 2 JL d (*2 

° l R l r 



FIGURE I -2(b) - BREAKDOWN OF THE TOTAL DRAG OF BLUNT AND SHARP 

AFTERBODIES INTO INTEGRATED VALUES OF SKIN FRICTION 
AND PRESSURE DRAG 



the use of the Direct Method is not recommended either for the purpose of data 
analysis or for the optimization study for the development of an airfoil con- 
figuration. 

The profile or total drag of the airfoil section is usually calculated by 
some form of the Indirect Method. The calculation of the profile drag of a 
given airfoil section by use of the Indirect Method is dependent upon the abili- 
ty of this method to reliably calculate the shape of the velocity profile at a 
very large distance downstream of the airfoil trailing edge. This very large 
downstream distance corresponds to the distance behind the trailing edge at 
which location static pressure is essential ly constant and > i ts value is equal to 
free stream static pressure. The drag of the given airfoil can then be calcu- 
lated by substracting the flow momentum corresponding to the above velocity 
profile at the large downstream distance from the corresponding momentum of the 
free stream air. 

Use of the Indirect Method for the calculation of the profile drag, however, 
requires knowledge of the characteristics of the boundary layer on both upper and 
lower surfaces of the given airfoil in the vicinity of its trailing edge. These 
boundary layer characteristics, known as initial conditions, can be calculated, 
for example, for attached flow conditions by the use of methods in Reference 1. 

The Indirect Method for calculations of profile drag, which is available up 
to the present time, is the Squire S Young's method, of Reference 2. The simpli- 
fied solution for the profile drag by the above method is obtained, however, at 
the sacrifice of making some assumptions which are inconsistent with physical 
reality. For example, logarithmic relation of velocity profile in wake, namely, 
the assumption of constancy of the ratio Log^U^/U^/H-l along the wake, U^jCon- 
tradicted by several experimental measurements in airfoil wakes. MoreoveryJ^me 
Squire £ Young approach can be used to advantage in some instances, its limita- 
tions prevent reliable application to many problems of practical interest. 

It is thus desirable to develop a generalized theoretical prediction method 
for the calculation of profile drag of an arbitrary airfoil, which is free from 
restrictions and limitations, such as discussed in the above paragraphs. Speci- 
fically, the method should be valid at high and low angles of attack and 
applicable to either single or mul ti -component airfoil. Specific objectives to 
be accomplished under the present studies are: 

(1) To obtain measurements of velocity profile and pressure distributions 
in the wake of a single component airfoil for the purpose of develop- 
ing a physical wake flow model. 

(2) To derive the equations to provide a generalized mathematical model of 
airfoil wake and to develop a numerical method and a computer program 
for the solution of these equations. 

(3) To determine the validity of the method by correlation with measure- 
ments of the viscous flow in the airfoil wake and profile drag 
measurements. 
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II. THEORETICAL STUDY 


11.1 General Discussion 


, The profile drag or total drag of the airfoil section is made up of two 
parts, namely, skin friction drag and pressure drag. Thus 

C[) profile ~ ^pressure + C °skin friction 

In the Direct Method, the total drag of the airfoil section is computed by 
evaluating the two components of the right-hand side of equation (I 1 - 1 ) separate- 
ly, namely, (a) calculation of pressure drag and (b) skin friction drag. A brief 
description of this method was given in Section I. In the Indirect Method, the 
total drag is computed without evaluating the two components of the right-hand 
side of equation (l 1-1). This is accomplished from the knowledge of boundary 
layer parameters on the upper and lower surface in the vicinity of the trailing 
edge. 


Figure 11-1 illustrates the basic principle behind the Indirect Method for 
the purpose of calculating the profile drag of single or mul ti -component airfoil 
section. From s the knowledge of the characteristics of the boundary layer on the 
upper and lower surfaces in the vicinity of the trailing edge, the Indirect 
Method is used for the purpose of calculating the velocity profile in the wake, 
at a large downstream distance behind the airfoil trailing edge, where static 
pressure distribution has stabilized to the freestream value. The expression 
for the profile drag of the two-dimensional airfoil section can then be written 
as 

+ oo 

Total Drag * sp j 11(1/^ - u) dy ( I I -2) 

— oo 

where s = span of airfoil 

u = local velocity in the wake boundary layer 
U = freestream velocity 
p = densi ty. 

The drag of the airfoil can then be determined if it is possible to use 
known velocity profile characteristics at the trailing edge of the airfoil to 
calculate the far downstream wake velocity distribution, y ) , point 4, 

Figure 11-1. 

U (x-*-», y) = F '( U (x = T.E., y)’ P l’ P 2’ P 3 •** P n) (N-3) 

where Pj , i =1 ... n are n different parameters required for the functional 
representation of equation (11-3). For example, in the case of a single compo- 
nent airfoil with a blunt trailing edge and with boundary layer separation ahead 
of the trailing edge in its vicinity, Pj may be free stream velocity, P 2 may be 
local outer flow velocity at separation, P 3 the trailing edge thickness and P 4 
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FIGURE II-1 - SCHEMATIC REPRESENTATION OF AIRFOIL WAKE FOR COMPUTATION OF PROFILE DRAG 




may be the distance between separation point and the trailing edge. Each pair 
of parameters are also related to each other through the general velocity dis- 
tribution U( x y) at any point in the wake. For example, if is the 

velocity at tfie edge of the boundary layer in the wake, and is the free 
stream velocity, then it is found from experimental data that it is possible to 
express the functional relationship as follows 

u e(x) - F (U.. D) . (M-*) 


where 


o 

Thus, (n-1) auxiliary equations, such as equation (11-4), are required for the 
complete solution of equation (I I — 3) - 

In the case of a single component airfoil with a sharp trailing edge, which 
does, not exhibit any separation ahead of the trailing edge, only two parameters 
in equation (11-3) need to be considered* These parameters are and the 
pressure distribution or U e , . along the wake. Squire and Young proposed the 
following empirical relat ion x 'for the purpose of obtaining an analytic expres- 
sion for the drag of a single component airfoil, i.e. 


o 

* 1 


u dy s average velocity in the wake boundary layer. 


H (X) - 1 


Log- (• 


e 'u 


(T.-.E-O. 


h T.E. ' 1 


Constant 


( 11 - 5 ) 


where 

Hfxl = wake boundary layer form factor 


U 


e (x) 


= velocity distribution at the outer edge of the wake. 


By making use of equation (1 1-5) in the von Karman momentum integral 
equation. Squire and Young were able to derive an analytic expression for the 
total drag coefficient of a single component airfoil with a sharp trailing edge. 
This expression, which is theoretically valid only in the absence of boundary 
layer separation ahead of the airfoil trailing edge, is given by 


U, 


Cd - 2(6 U + e L ) (- 


«T. E. 


(H t . e . + 5)/2 


u 


•) 


( 11 - 6 ) 


It may be pointed out here that the above simple analytic expression for 
drag coefficient is obtained at the sacrifice of making some simplifying assump- 
tions which are inconsistent with physical reality. For example, experimental 
data in the wake behind airfoils indicate that the empirical relation of equa- 
tion (ll~5) used by Squire and Young does not hold. Figure 11-2 shows the plot 
of experimental data for the ratio (Log e (U /U e ))/(H-1) at various distances aft 
of the trailing edge for a Joukowski airfoTl. These results show that 
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assumption Implied by equation ( 1 I — 5) in the Squire and Young method is un- 
realistic. Specifically, Squire and Young's approach is unable to deal 
directly with the following situations: 

(1) Blunt trailing edge airfoils. 

(2) Airfoils composed of two or more components. 

(3) High or low angles of attack when boundary layer thickness on the 
upper and lower surfaces of the trailing edges are very difficult. 

(A) Incipient boundary layer separation in the neighborhood of the 
trailing edge on either the upper or lower surface or both. 


11.2 Description of the Present Theoretical Method 


Figure I I “3 shows the curvilinear system of coordinates in which x denotes 
the distance along the locus of minimum velocity and y the distance normal to it. 
The radius of curvature of the lower edge of the wake will be denoted by Ri and 
of an equidistahce line by R. It is assumed in this theoretical development that 
Rj is large in comparison with the width of wake, but not very large in compari- 
son with typical wake layer width. Defining curvature K as 


K = 


Ri 


then one has 


R 


1 

1 ± Ky 


(1 1-7) 


( 11 - 8 ) 


The complete Navier-Stokes equations for the orthogonal directions, when 
local curvature in flow is taken into account, can be written as follows: 


jL 3u 3u uv 

R U 3x V 8y R p R 9x 


= V 


afu 

r 2 ax 2 


a 2 u 

ay 2 


R i 3 8 ( rV au 1 3u 
r 3 y 3x * 3x R 3y 


R2 R 2 <* x R 3 9x 

T 

X - Momentum Equation 


(11-9) 
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( 11 - 10 ) 



3v . 3v 

u — — + V 

3x 3y 


“i + llL = v R l 2 8 2 v , 3 2 v 

R P 3y ' R 2 3x 2 3y 2 


_1_ J_ , J_x 3v 1 3v 

r 3 Y 3x 'R ' 3x R 3y 


v 

R 3 



_a_ 

ax 



2 


^1_ 3u_ 

R2 3x 


Y Momentum Equation 


The equation of continuity, considering effects of curvature, can be 
written as 


3u + 

ay ay 




= 0 


( 11 - 11 ) 


By performing time averaging operation on the terms of equation (I 1-9) and 
(l 1-10), the equations for the mean turbulent flow In the wake behind the trail- 
ing edge of the airfoil can be derived. These equations can be further simpli- 
fied from the consideration of the order of magnitude analysis. The following 
time averaged equations can be written after appropriate simplification for the 
presently considered wake flow. 


u 


ao r 

ax Ri 



1 9P + R_1 3l££_ + _Ll Txy 

P 3x Rj p 8y Ri P y 


( 11 - 12 ) 


15P = u£ / 1 \ _ _R_ 3Txy 

p 3y R M + y/R' R 3x 


(11-13) 


£ + = 0 • 

In the above equations, terms of the order 6/R and of higher magnitude are 
retained and the terms of the order smaller than 6/R^ are eliminated; here 6 is 
the characteristic width of the wake boundary layer. The shear stress T(x,y) 
which appear in above equations is the algebraic sum of laminar and turbulent 
shear stress contributions, i.e. 


T(x,y) = ]i |y- - p vC v 1 . (11-15) 

In the further derivation of equations for the wake flow the bar above 
qualities, such as u, v, p, etc., will be omitted; however, it will be implied 
that they represent the time averaged value at any point in the flow behind 
airfoil wake. 

Physical Flow Model : Figure I 1-4 shows the physical flow model of the wake 

aft of blunt base mult? -component airfoil. This flow model was derived from 
available experimental data. The success of predicting the complex viscous flow, 
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FIGURE II-4 - PHYSICAL MODEL FOR THE VISCOUS FLOW IN THE VAEE OF SINGLE COMPONENT OR WO-COMPONENT AIRFOILS 


















(Continued) 



such as shown in:Figure 11-4, for an arbitrary airfoil case by direct solution 
of simplified Navier-Stokes equation (11-12), (11-13) and (11-14) would be 
highly improbable, if not impossible. For this reason, the flow model of 
Figure 11-4 is divided into various regions and layers in which the flow be- 
havior has certain physical characteristics. For example, experimental 
measurements in the channel flow suggest that Region 1 can be divided into 
three Zones LI, L2, and L3. In Zone LI, the maximum circulatory flow is 
present as the pressure is almost constant, whereas in Zone L2, pressure in- 
creases rapidly to near stagnation value at the point of confluence. In Zone 
L3, which is infinitesimally small, pressure decreases at a rapid rate and the 
changes in velocity profile shape are quite abrupt. 

From the present experimental measurement in the wake of a sharp trailing 
edge airfoil, it is observed that velocity profiles in the individual layers 
of Region IV of Figure 11-4 become "similar 11 or "one parameter family" if 
appropriate parameters are chosen. If the characteristics of the flow are 
such that it is possible to derive similarity parameters for all layers in 
different regions of Figure 11-4, then sets of ordinary differential equations 
for this wake flow can be derived from the governing partial differential equa- 
tion such as (11-12), (M-13), and (ll-l4). These sets of ordinary differential 
equations, however complex, can then be solved by presently available numerical 
methods for the solution of the flow in the wake of an arbitrary airfoil. The 
discussion of similarity parameters for velocity profiles and generalized 
physical parameters for shear stress and pressure distributions in the wake be- 
hind single component airfoil will be given in Section IV. However, in the 
following derivations of theoretical equations for wake flow behind thin trail- 
ing edge single component airfoil section, it will be assumed that such simi- 
larity or generalized physical parameters do exist for the considered viscous 
flow. 


11.3 Theoretical Derivations 

Equations for Zone L3 of Region 1 : As the presently considered case is for 

single component airfoil with the sharp trailing edge, Zones LI and L2 are 
absent; solution for the flow is thus required for Zone L3 only of Region 1 for 
this case. Schematic representation for Zone L3 for single component airfoil is 
shown in Figure 11-5- Initial conditions for the boundary layer quantities on 
upper and lower surfaces of airfoil are specified at Section AA, shown in this 
figure; this Section AA corresponds to either trailing edge location or is 
situated at the point of incipient separation on the airfoil which is in the 
neighborhood of airfoil trailing edge. The Section BB, shown in Figure 11-5, 
corresponds to the end of Zone L3 of Region I. This Section BB is located at an 
infinitesimal distance from Section AA; however, due to rapid mixing of boundary 
layers of upper and lower surfaces of an airfoil in Zone L3, characteristic 
viscous flow quantities such as static pressure, displacement and momentum 
thickness and magnitude of minimum velocity U w change at much faster rates in 
Zone L3 than in downstream regions of flow. A starting value of is obtained 
by making a momentum balance at the trailing edge. The use of momentum 
integral equation in the following form is made: 
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Station A-A — Beginning of Zone L-3 
Station B— B — End of Zone L-3 
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FIGURE II-5 - SCHEMATIC REPRESENTATION FOR ZONE L-3 FOR SINGLE COMPONENT AIRFOIL 


7P+ (H + 2) ~ = 0. 

dx U e dx 


( 1 1 - 16 ) 


The parameter and function for the similarity of velocity profiles in Zone L3 
are assumed to have the following form, namely, 


For Upper Half Wake: 


y - Y2 


fO^) “ M 


U e .. - U 


V3-Y2’ U e u - % 


(11-17) 


For Lower Half Wake: 


y 2 ~ y ,, , u 

n 2 y 2 - Yj* f(n 2 } Ue L - U w 


( 11 - 18 ) 


By making use of equations (11-17) and (11-18) in the usual definitions for 
momentum and displacement thickness, the expressions for these quantities for 
upper and lower wakes in Zone L3 can be derived in terms of parameters of equa- 
tions (11-17) and (11-18). Thus, 


eUpper 

Wake 


0 Lower 
Wake 


r u Wi f 1 

" (y 3“ y2) huT" 0 ] f(n i )dn i 
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e u 
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fln^dnl 

0 J 


(11-19) 


(y 2 -yi) - 1) f(n 2 )dn 2 + (^ - i) 2 | f(n 2 ) 2 dn 2 j ( 11 - 20 ) 


Si S r =<Vs-v 2 ) 0-A ff(Vdn 


(11-21) 


<5*Lower 


= (y 2 -yj) 0 - ir 1 ^ I 


2 2 


( 11 - 22 ) 


By making use of equations (11-17) through (11-22) in equation (11-16) and 
performing needed mathematical manipulations, the following quadratic equation 
can be derived for computing the ratio (U^/Uga) at the end of Zone L3. 

U Wl 2 Uu 

A 1 (jj-^) + A 2 Cjj-4 + A 3 = 0. ( I 1-2: 
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Where the vaJues of coefficients are determined from following equations. 
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The variation of the width of the upper and lower wakes in Zone L3 is ex* 
pressed by the following growth rate equations. 


and, 


z!r <Y3-y 2 ) = ci u • 



■h (y 2 -yi) = ci L . 



(11-24) 


where Cl u and C 1 are constants and initial values of (y 3 “ Y 2 ) and (Y2 “Yi) f° r 
equations (11-24) and (11-25) correspond to values of boundary layer thickness at 
Section AA on upper and lower surface of the airfoil, respectively. The implied 
assumptions for the validity of equations (11-24) and (11-25) are that growth 
rate of these layers is controlled by the transverse perturbation velocity and 
that the transverse perturbation velocity is proportional to the average gradi- 
ents of velocity for upper and lower wakes in Zone L3* Simultaneous solutions 
of equations (11-23), (11-24), and (11-25) give the starter values for (Uy /U ei| ), 
(y3~y2) and (Y 2 ~Yi) at the end of Zone L3 of Region 1 . In the case of 
single component airfoil, the Regions II and I I I of Figure 11-4 are not present; 
and, thus, these starter values form the initial conditions for Region IV in the 
wake behind single component airfoil. 

Equations for Region IV : For single component airfoil, the Region IV ex- 

tends from the end of Zone L3 of Region I for a distance very far downstream of 
the airfoil trailing edge where the pressure has stabilized to the free stream 
value. Experimental data of the velocity profiles for layers 8-5 and 1-8* in 
Region IV indicate that velocity profiles in these layers become ''similar 11 if 
the similarity parameters and similarity functions are defined in the following 
manner: 


For layer 1-8: 


Y8 “ Y 

n3 '9T r T7 ; f(n3> 


Ue L - u 

u e L " U W8 


where Ul = velocity at any point y in layer 1-8 

UeL = velocity at lower edge of wake boundary layer 

Uy 8 = minimum velocity on the locus of Ys( x ) 

y £ = distance y in the layer 1-8 where u -1/2 (Uy 8 + U 6 |_) 


( 11 - 26 ) 
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(M-27) 


and for layer 8-5 'define: 
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y~y 8 
y 8 c “ y s ; 


f (114) 


Ue u - u 

Ue u - Uwa 


where u *» velocity at any point y In layer 8-5 

Ue u “ velocity at upper edge of wake boundary layer 

y 8 g ■ distance y in the layer 8-5 where u ■ 1/2 (Ue u + Uy 8 ). 

In the definition of s i mi lari ty parameters r> 3 and 04. as given by equations 
( 1 1—26) and (ll-27)» the distances y X Q and y 0 Q are used because it is difficult 
to determine the exact values of boundary layer thickness from measurements at 
large distances downstream of the airfoil trailing edge. The Euler equation at 
the upper edge of the wake can be written as: 


dp d u e u 

dx p e u dx 


(11-28) 


whereas for the lower edge of wake, the Euler equation is written as: 


dP _ .. d U e|_ 

d7 = ‘ Ue L ~dT 


(I 1-29) 


The expression for the transverse velocity V(y) at any point in the layer 
1-8 or layer 8-5 can be derived from the continuity equation (11-11) as follows, 


V(y) S& 1 - v<y«> Hr - J 


(11-30) 


y8 


where R(y) » radius of curvature of stream line in wake boundary layer at any 
point y 

R2 * radius of curvature for the locus of minimum velocity y = Y8(x) 
R1 ■ radius pf curvature for lower edge of wake boundary layer 
V(y8) - transverse velocity on the locus y = yg( x ) of minimum velocity 

U 8(x)* 

Applicable boundary layer conditions for the layer 8-5 and layer 1-8 are: 
at y « y 5 : u » Ue u , f(n 4 ) - n 4 - K lu > Cp * Cp(y 5 ), t * x(y 5 ) * 0, 


and — 

P 


1 d Pe u 


- u e 


d Ue., 


dx Je u ~dT 

at y - y e : u - Ue 8 , f(n4) - f(n 3 ) - 1 , 113 

and V(y) - V(y 8 ) 
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0, Cp - Cp (y 8 ) , 


23 


and at y = y u = U eL , f(n 3 ) = 0 , n 3 = K lL 

Cp = Cp(y x ) , t = x( Yl ) = 0, and 1 = - UgL (11-31) 

In order to obtain momentum integral equation for layer 8-5, individual 
terms of equation can be integrated from y = y 8 to y = y 5 . During the mathe- 
matical manipulat ions, use is made of equations (M“ 30 ), ( 11 - 29 ), ( 11 - 26 ), 
boundary conditions (11*31) and Leibnitz's rule. Radius of curvature R for the 
wake flow boundary is assumed large compared to the thickness of boundary layer 
for the purpose of analytical simplification. Integration of the term 

y r 5 

' sr dy 

y 8 , 


needs some mention at this time. I h order to perform this integration, the 
information regarding variation of integrand in the layer 8-5 is necessary. 
Experimental data indicate similarity in ve loci ty prof i le for longitudinal 
velocity u; however, such general i zed information regarding transverse compo- 
nent v is not avai lable-at the present time. Under this c i rcumstance , -the 
results of analytical solution for Tollmien's (Reference 3) plane turbulent 
source are used for the purpose of evaluating this integral. Figure ll- 6 (a) 
shows the schematic representat ion of Tollmien's plane turbulent source. 

Tollmien investigated the velocity distribution in the free jet flow which is 
formed when the flow of fluid emanates through the narrow rectangular slit. The 
x-axis coincides with the axis of symmetry and the y-axis is perpendicular to 
x-axis at the exit of the slit. t The results of analytical solution, obtained by 
Tollmien for the flow conditions of Figure I I - 6 (a ) are used for the present 
purpose as follows: 
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wnen y = y 8 , 


Substitution of (11-33) in (11-32) and simplifying, 
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Figures 11-6 (b) and I I -7 show the plots g 2 (?) and gi (£). g 2 (C) versus (?) ; 
these curves are obtained from theoretical results for Toll mien's plane turbu- 
lent source solution. From the above figures, the values of integrals can be 

written as: 

K 2U K 2u 

| g 2 (S) dC : - -0214 and J gi (5) g 2 U) d£ = .0065 (I 1-35) 

O 0 


Thus, the momentum integral equation for the layer 8-5 can be derived as 
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(11-36) 


If we assume that the rate of growth of the layer 8-5 is controlled by the 
transverse perturbation velocity, and — as experimental evidence indicates — 
similarity of velocity profile for this layer, then this growth rate can be 
expressed by the following equation: 
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where C 2u is an empirical constant. 

In the manner similar to the derivation of equation (11-36), the momentum 
integral equation for the layer 1-8 can be derived as: 
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R 8 _5 and Rj_ 8 , appearing in equations (11-36) and (11-38) can be approxi- 
mated as the arithmetic mean values of radii of curvature for the layers 8-5 and 
1-8 respectively. The growth rate equation for the layer 1-8 can be written in 
a manner similar to that for layer 8-5> as follows: 
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where C 2 ^ is an empirical constant. 

Thus, equations (11-36), (1 1-37) » (11-38) and (11-39) are four equations 
for the simultaneous solution of four variables yj, y 2 , and Uu g in the Region 

IV, which is shown in Figure 11-4. These equations can be arranged in a manner 
to form the Initial Value Problem. Various single-step or multi-step methods 
can be used for the solution of this initial value problem. The most commonly 
used methods are the single step Euler method, the single-step modified Euler 
method, the multi-step predictor corrector method, the Runge-Kutta method, and 
a few others. The choice of a method depends upon the particular problem and 
is governed by the desired accuracy, time of computation, core size available in 
a particular computer, and other factors. 


Equations (11-23), (11-36), and (11-38) contain terms such as 


Kl 

f (n)dn , 

J 

0 


Ki 

f f 2 (n)dn, 


K 2 

» 

9l (C)g 2 (5)d5» 

J 

O 


shear stress terms, terms containing shear and pressure integrals and CP dis- 
tributions at the edges of various layers. In order to be able to solve these 
equations, however, the values of the above quantities either have to be known 
priori and/or more auxiliary equations are required which express the above 
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quantities as the functional relationships in terms of dependent variables. As 
the viscous flow under consideration is turbulent wake flow, theoretical expres- 
sions for the above quantities are not available as in the case of laminar 
boundary layers. Recourse is then made to experimental measurements to obtain 
empirical expressions for the above parameters by the use of experimental data 
for the particular flow which is being investigated. This matter is further 
discussed in Sections III and IV. 


III. EXPERIMENTAL WORK 


In order to facilitate understanding of turbulent wake flow, which exists 
behind the trailing edge of airfoil model, for the purpose of developing an 
analytical model and also to check the validity of theoretical predictions, an 
experimental program was conducted in the research wind tunnel facility at 
Lockheed-Georgia Company. In this section, a brief description is given for 
airfoil model, the experimental facility, probes for measurements of velocity 
and static pressure profiles, and types of tests. Detailed description of 
instrumentation, sidewall blowing requirements in tunnel working section, and 
data acquisition and reduction is given in the appendices. 


I 11.1 Description of Airfoil Model 

Airfoil model is 15 percent thick and is symmetrical airfoil section. It 
is basically the NASA four digit airfoil section with slight modification near 
the leading edge. Figure I I 1-1 shows the geometry of the test airfoil configu- 
ration; the difference in geometry near the leading edge of the test airfoil 
from NACA 0015 is also indicated in this figure. Computer program of Reference 
1 indicated the existence of laminar stall on the original NACA 0015 airfoil at 
a Reynolds number of approximately 1 million. The leading edge geometry for 
the test airfoil was arrived at by the use of the computer program of Reference 
1 such that turbulent separation near the trailing edge was predicted before 
the occurrence of laminar stall prediction on this test airfoil. The above 
phenomena then makes it possible to conduct systematic studies for the flow in 
the wake behind the airfoil trailing edge. 

Figure I I 1-2 shows the comparison between predicted C[__ a curve of the test 
airfoil with experimental measurements. It can be seen from this figure that 
viscous prediction agrees quite well with experimental measurements up to an 
angle of attack of a * 9° to 10°. Experimental measurements of boundary layer 
velocity profiles indicate that turbulent boundary layer separation on the 
upper surface of the test airfoil near the trailing edge appears for angles of 
attack a greater than 9° to 10°. Theoretical methods which are used in Refer- 
ence 1 for computations of aerodynamic characteristics of airfoils are not valid 
in the presence of separated flow on the airfoil surfaces. For this reason, 
theoretical predictions of lift coefficients differ from experimental measure- 
ments for angles of attack a greater than approximately 9° to 10°. 


31 









Pressure distribution for the test airfoil at various angles of attack 
are shown in Figures lll-3(a)» lll~3(b) f III —3 (c) and lll-3(d). Computed 
results for potential and viscous pressure distributions are shown in the above 
figures and experimentally measured pressures at various chordwise locations on 
test airfoils at various angles of attack are shown plotted for comparison 
purposes. Comparison shown in Figures ! I I -3 (a) and lll-3(b) indicate that 
theoretical predictions agree quite well with experimental measurements for 
angles of attack a equal to 0° and 8°; however, for angles of attack a of 12° 
and 14°, agreement between computed and experimental pressure distributions is 
not so good because of the existence of turbulent boundary layer separation on 
the upper surface of the airfoil. 

Figures lll-4(a) and lll-4(b) show the plots of computed boundary layer 
development on the upper surface of the test airfoil for angles of attack of 8° 
and 12°, respectively; experimental measurements of boundary layer quantities on 
the airfoil surface aft of X/C of 0.8 are also shown in these figures. For the 
case of angle of attack of 8°, shown in Figure lll-4(a), incipient separation 
may be present aft of X/C of O .96 near the trailing edge of the airfoil (as 
indicated in this figure by the computed value of the turbulent boundary layer 
form factor H which becomes greater than 2.0 near the trailing edge). In this 
case the computed values of boundary layer displacement and momentum thicknesses 
and form factor agree reasonably well with experimental data. However, when the 
angle of attack is 12°, the results of computations of boundary layer quantities 
do not agree very well with experimental measurements on the upper surface of 
the airfoil because of the presence of appreciable region of flow separation on 
it. The boundary layer separation takes place at X/C £ 0.8 at an angle of 
attack of 12° approximately, and for this reason measured boundary layer thick- 
nesses are much larger than computed values aft of X/C = 0.8. 

The model was mounted vertically in the working section of the tunnel 
spanning the 76.2 centimeters test section. The chord of the airfoil model is 
29.2 centimeters. Total of 16 static pressure orifices are provided on the 
upper surface and 10 orifices are provided on the lower surface of the model 
along Butt Line 0.00, i.e. along the center span of the model. In addition, a 
total of 15 additional static pressure orifices are provided on the upper and 

lower surfaces at Butt Lines ±3.00 for checking the two dimensionality of the 

flow on the airfoil model. Table 3” 1 gives coordinates and Table 3“2 shows 
the chordwise locations of the static pressure orifices on the model. 

Forward of chordwise location X/C of 0.7> the model is comprised of a 
central load-carry ing beam covered by a contoured maple shell. The shell con- 
sists of a solid maple leading edge and upper and lower surface ‘‘skins 11 . The 

leading edge is bolted to the front of the beam. Aft of X/C of 0.7 the model 

is made of solid aluminum. A tolerance of ±.0127 centimeters is maintained on 
the airfoil contour across the span. 


III. 2 Wind Tunnel Facility 

Experiments for the studies of the wake flow behind the single component 
test airfoil were conducted in the Lockheed Research Wind Tunnel facility which 
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FIGURE III-3(a) - COMPARISON OF PRESSURE DISTRIBUTIONS FOR 

THE TEST AIRFOIL (SHARP T.E. ) AT AN 
ANGLE OF ATTACK OF 0° 
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FIGURE III-4(b) - BOUNDARY LAYER DEVELOPMENT ON THE UPPER SURFACE OF 

THE TEST AIRFOIL (SHARP T.E.) AT AN ANGLE OF 
ATTACK OF 12* AND COMPARISON WITH EXPERIMENTS 




Table 3=1 - COORDINATES OF THE TEST AIRFOIL 
SECTION (SYMMETRICAL AIRFOIL) 


Airfoil Chords, 

Stations, 

Percent Chord 
X/C x 100 

0.00 

0.23 

. °°57 ,,, . 

1.14 

1.72 
2.29 
5.43 
4.58 

5.72 
6.87 
8.01 

10.30 

12.59 
14.88 

17.17 

20.60 

25. 18 
30.00 

- 34.33 

38.91 
43.49 
48.07 
52.65 
57.22 

62.95 

68.67 

74.39 

80.11 

85.84 

91.56 

96. 14 

100.00 


= 29.2 inches 

Airfoil Section Ordinate 
Percent Chord 

Z/C x 100 

0.00 

1.45 
2.20 
2.87 
3.36 
3.76 
4.33 

4.73 

5.00 

5.32 

5.45 
5.90 

6.32 
6.64 
6.93 
7.21 

7.46 
7.50 
7.44 
7.29 
7.04 
6,72 

6.33 
5.90 

5.38 

4.78 

4.08 

3.35 

2.52 

1.60 

0.77 

0.00 



Table 3-2 - LOCATION OP PRESSURE ORIFICES 
ON THE AIRFOIL MODEL SURFACES 


Tube 


Number 

Surface 

Butt Line 

Z/Q 

1 

Upper 

0.0 

0.00 

2 

Upper 

0.0 

0.01 

3 

Upper 

0.0 

0.05 

4 

Upper 

0.0 

0.15 

5 

Upper 

0.0 

0.40 

6 

Upper 

0.0 

0.60 

7 

Upper 

0.0 

0.65 

8 

Upper 

0.0 

0.69 

9 

Upper 

0.0 

0.74 

10 

Upper 

0.0 

0.78 

11 

Upper 

0.0 

0.82 

12 

Upper 

0.0 

0.86 

13 

Upper 

0.0 

0.90 

14 

Upper 

0.0 

0.93 

15 

Upper 

0.0 

0.96 

16 

Upper 

0.0 

0.985 

17 

Upper 

+3.0 

0.05 

18 

Upper 

+3.0 

0.40 

19 

Upper 

+3.0 

0.69 

20 

Upper 

+3.0 

0.93 

21 

Upper 

+3.0 

0.985 


imber 

Surface 

Butt Line 

x/c 

22 

Upper 

-3.0 

0.05 

23 

Upper 

-3.0 

0.40 

24 

Upper 

-3.0 

O .69 

25 

Upper 

-3.0 

0.90 

26 

Upper 

-3.0 

0.96 

27 

Lower 

0.0 

0.02 

28 

Lower 

0.0 

0.04 

29 

Lower 

0.0 

0.10 

30 

Lower 

0.0 

0.35 

31 

Lower 

0.0 

0.50 

32 

Lower 

0.0 

0.74 

33 

Lower 

0.0 

0.80 

34 

Lower 

0.0 

0.90 

35 

Lower 

0.0 

0.95 

36 

Lower 

0.0 

0.985 

37 

Lower 

+3.0 

0.35 

38 

Lower 

+3.0 

0.985 

39 

Lower 

-3.0 

0.35 

40 

Lower 

-3.0 

0.90 

41 

Lower 

-3.0 

0.985 
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is located in the Aerospace Sciences Laboratory of the Lockheed-Georgia Company, 
The cross-section of the test section is rectangular having an area of O.836 
square meters with a height to width ratio of 0.7; this results in a test sec- 
tion of approximately 76.2 x 109.2 centimeters with a length of 122 centimeters. 
This wind tunnel is a closed circuit, single return, low speed wind tunnel and 
.is powered by a 400 horsepower synchronous speed induction motor driving a 1 83 
centimeters, 10-bladed, single stage axial flow fan. The fan speed can be varied 
over .the range from 0 to 1150 revolutions per minute by means of an eddy current 
clutch. A standard friction brake rated at 103.7 meters kilogram torque will 
bring the fan to rest from full speed in less than 15 seconds. 

The velocity range of the tunnel with an empty test section is 0 to 91*5 
meters per second resulting in a maximum dynamic pressure of 508.3 kilograms 
per square meter. The dynamic pressure variations in the test section can be 
maintained within approximately 0.1 percent during a boundary layer survey 
which typically takes about ten minutes in this test facility. With the present 
screens, the turbulence factor is less than 1.3 based on turbulence sphere 
measurements. At a dynamic pressure of 293 kilograms per square meter, the 
nominal Reynolds number is .439 x 10 6 per meter of chord length. A general lay- 
out of the test facility is shown in Figure 1 1 1 — 5 * 

Figure 1 1 1-6 shows the variation of static pressure along the test section 
centerline at various free stream dynamic pressures. This figure shows that 
the static pressure at the location of the quarter' chord of the model is within 
3.bout 3 percent of the static pressure at any other location on the model center- 
line. The calibration done at the location of the model quarter chord was used 
In reducing data from these tests. This calibration is shown in Figure ill-7. 

The test section is equipped with blowing boundary layer control slits on 
the ceiling floor. The auxiliary air, for boundary layer control on the test 
section side walls, is provided from the storage tanks which are capable of 
delivering air at a maximum pressure of 22.73 kilograms per square centimeter. 

A regulator reduces this pressure to 11.23 kilograms per square centimeter, a 
few meters from the test section. This air is further throttled to a desired 
pressure in the plenums for boundary layer control slits. 

The boundary layer control slits, which are located on the ceiling and the 
floor of the test section, are 77.2 centimeters wide and .0635 centimeters high. 
These slits are located approximately one chord length upstream of the model 
leading edge. The heights of the blowing slits can be varied from .0254 to .254 
centimeters. Relatively high velocity air is introduced through the slits into 
the tunnel floor and ceiling boundary layers such that boundary layer separation 
is prevented on these side walls when the given model. is in the tunnel at a 
given angle of attack with respect to the free stream. When the angles of 
attack for a given two-dimensional model in the tunnel is increased from low 
values to high values corresponding to condition, the boundary layer sepa- 

ration on the tunnel side walls occur upstream of the model leading edge due to 
the presence of adverse pressure gradient in the flow direction. If the above 
flow separation is present, then the reliable two-dimensional data cannot be 
obtained, especially near the Cu,** conditions. The method used to determine 
the proper amount of blowing and the desired pressure ratio required to elimi- 
nate this boundary layer separation on wind tunnel side walls is discussed in 



FIGURE III- 5 - TEST FACILITY GENERAL LAYOUT 
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Appendix A. Figure I! 1-8 shows the details of the boundary layer control 
system used for the present study contract. 


I 1 1-3 Special Instrumentation - Pressure Probe, Hot Wire 
Anemometer and Probe Drive Mechanism 

Measurements of flow in the wake as well as on the upper surface of the 
test airfoil in the vicinity of the trailing edge were performed by the use of 
pressure probe and hot wire anemometer probe. In the case of sharp trailing 
edge airfoil, flow measurements in the airfoil wake were performed by use of a 
special type of pressure probe, whereas in the case of test airfoil with 1 
percent thick trailing edge, hot wire anemometer probe was used for wake flow 
measurements. In the case of blunt trailing edge, the velocity profiles have 
reverse flow or negative velocity associated with them up to some distance 
downstream of the trailing edge. Hence, in this case, pressure probe alone can- 
not be used for determination of flow quantities in this region and hence addi- 
tional instrumentation, such as hot wire probe is necessary for the determination 
of flow associated with reverse flow velocity profiles. The following paragraphs 
give the brief description of these two types of probes used during the experi- 
mental program of the present studies. The details of the additional instrumen- 
tation, used for acquisition and reduction of experimental measurements is given 
in Appendix B. 

1 1 1.3.1 Pressure Probe : Measurements of static and total pressures in the 

boundary layer, for the purpose of determining boundary layer velocity profiles in 
the wake of test airfoil as well as on its upper surface, were made using the 
specially designed three element probe shown in Figure 1 1 1 — 9 - For the purpose of 
developing accurate mathematical model for the calculation of wake flow and 
profile drag of airfoils, it was desired to know the true variations of static 
pressure distribution across the boundary layer of the airfoil wake at various 
chordwise locations. Hence the selection of the probe tip design for the 
measurements of total and static pressures in the viscous flow in the wake was 
based on the following considerations. 

It is known that it is relatively easy to measure an accurate value of total 
pressure in the flow with li.ttle effect arising from incidental variation in flow 
conditions. On the other hand, the selection of hole positions on the probe for 
the measurements of static pressure in the flow is much more difficult. This is 
so because it is observed that the measured value of the static pressure in the 
flow is very sensitive to any parameters, e.g. pitch angle of the probe, local 
flow direction in the flow, conditions of the boundary layer on the probe, prox- 
imity of static pressure holes on the probe to solid boundaries and several other 
factors. The above difficulty can, however, be overcome if it is possible to 
derive the value of velocity in the boundary layer from measured pressure dif- 
ferences arising between holes which sense pressures greater than the local 
static pressure but which show negligible sensitivity to the above-mentioned 
parameters. A simple relationship for dynamic pressure then still holds good and 
can be expressed as. 
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FIGURE III-9 - TOTAL-STATIC PRESSURE PROBE COMBINATION 
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where 


j pu 2 = K (P x - P 2 ) 


(MI-1) 


P 


u 

?l and P 2 
K 


dens I ty 

local velocity 

pressure sensed at different hole positions, and 
calibration factor and may depend on wind speed. 


Probe tip design which has been fabricated in order to remove the above- 
mentioned sensitivities, and hence to be able to measure local velocity and 
static pressure to desired accuracy are shown in Figures MI-9 and 111-10. The 
upper and lower right-hand elements are designed to measure both static and 
total pressures. The lower left-hand element is designed to measure total 
pressure only. The total pressure tubes on the lower element thus provide a 
check or comparison. The upper element is used when the thickness of boundary 
layer exceed 7-62 centimeters. 

The left-hand, or total pressure, element is constructed of 0.0508 centi- 
meter diameter stainless steel tubing which has been flattenee at the end as 
shown in Figure 111-10. Measurement of total pressure within .01015 centimeters 
of the model surface are possible with this arrangement. The two side orifices 
were tied together and used to measure static pressure in the boundary layer. 

They were formed by beveling the side tubes 30 degrees as shown in Figure 111-10. 


The beveled tubes on the upper and lower elements measure the pressure 
which is proportional to the difference between total and static pressure and 
this measured pressure is related to local dynamic pressure by a calibration 
factor. Thus, a calibration of the probe reading versus true static and total 
pressure is necessary. Calibration of the probe was accomplished in the wind 
tunnel using the pitot tube shown in Figure 111-11 as a reference. The cali- 
bration curves for this probe are shown in Figure 111-12. The above calibration 
was performed in the absence of the solid boundary and therefore the correction 
due to the presence of wall proximity has to be used when velocity profile 
measurements near the solid surface are desired. It has been found by MacMillan 
(Reference h) that, for a pitot tube of ratio d/D * 0.6, the measured value of 
total pressure is less than the true value when the probe is nearer than 2D to a 
solid boundary; here d is equal to internal diameter of the tubing and D is 
equal to outside diameter of the tubing for the total pressure tube. This wall 
proximity correction curve, which is obtained from Reference h , is shown in 
Figure 111-13- Thus, total pressure, static pressure, and dynamic pressure at 
a desired point in the boundary layer is determined from measurements by probe, 
which is shown in Figure I 1 I - 9 , as given by the following expressions: 


P T C ■ P T ^ P T, 


wc 


= 7 pu 2 - K Ci (P T ^ - P 12 ) 

p 12) 
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= 4.88 KILOGRAM/SQ. METER 



FIGURE III— 12 - BOUNDARY LAYER PROBE CALIBRATION 





where 


and 


P T = corrected value of total pressure 
‘c 

Py = total pressure measured by total pressure tube P and 
corrected by the use of Figure 111-12 

Py = wall proximity correction from Figure 111-13 
wc 

P$ c * corrected value of static pressure 

K = calibration factor for obtaining static pressure from measured 
pressure by beveled tube P 12 by the use of curves in Figure 

I I 1-12 

u = local velocity in boundary layer 

R = universal gas constant 

T = absolute temperature in boundary layer. 


111.3*2 Hot Wire Anemometer Probe : The hot wire anemometer is a very desir 

able instrument for the analysis of the micro structure of the velocity of a 
streaming gas or liquid. Measurements of the micro structure of a stream require 
a very small sensitive element having a short response time, sufficient sensi- 
tivity and little disturbing effect on the original stream. The sensitive 
element is a thin electrically heated wire suspended between two needle points 
or a thin electrically heated metal film fused to a glass support. Figure 
111-14 shows the sketch of hot wire anemometer probe; also shown in this figure 
is the enlarged view of the single sensor and sensor support used during the 
present i nvest i gat ion. 


The use of a hot sensing element for measurements of particle velocity in 
fluid flows relies on laws governing convective heat transfer. These laws are 
generally too complicated to permit a theoretical calculation of the relation 
between the particle velocity and the heat flux from a probe, and the relation 
must therefore be found experimentally using laws of similarity. The expression 
for the heat flux for the two-dimensional heat transfer in an incompressible 
potential flow is given by 


Q = K* {1 + / (27rpC p dU/K)} (T-T 0 ) 


where 


Q 

K 

£ 


§ 

u 

T 

To 


hea-t flux 

thermal conductivity 
length of wi re 
densi ty of fluid 

constant pressure specific heat of fluid 

diameter of wire 

flow velocity 

film temperature 

fluid temperature. 


(M 1-3) 


For thermal-equilibrium conditions, the rate of heat loss from the hot wire 
must be equal to the heating power generated in the wire by the electric current, 
i.e. it must equal to l 2 R/J. Thus, from the point of view of hot wire anemometry 
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FIGURE III - 14 - HOT WIRE ANEMOMETER PROBE WITH SINGLE SENSOR 




the relation between fluid particle velocity and electrically generated heating 
power takes the following form: 


R 

R- Rq 


I 2 


A + B /U 


(MI-4) 


where A and B are constant for a specific probe wire operating in a specific 
fluid. 


u = local velocity 
I = heating current 
R = operating resistance of hot wire 
Ro - resistance of hot wire at fluid temperature. 


However, empirical formula for the heat loss from a hot wire is generally 
written as, 


I 2 = A + B (u)° (MI-5) 

and the values of A, B and n are determined from calibration. The output of the 
anemometer is a bridge voltage V, and the squared voltages V 2 and V Q 2 (V Q = 
bridge voltage when the flow velocity is zero) are linearly related to the heat 
loss of the wire at the velocity in question and at zero velocity respectively. 

At very low velocities, the heat loss due to forced convection is a minor part 
of the total heat loss, resulting in a typical calibration curve shown in Figure 
~MI-15(a). In cases where measurements of high degrees of turbulence are con- 
cerned, such as in wake flow behind airfoil, or where large series of measure- 
ments are to be performed over a wide velocity range at a great number of points, 
such as in the case of boundary layer measurements over the surface of an air- 
foil, it is necessary to linearize the anemometer output voltage in order to 
avoid nonlinear distortion and ensure rational operation. Figure 111—15 (b) 
shows the plot of anemometer output voltage versus flow velocity made in con- 
junction with hot wire anemometer. 

There are two modes of operation for the hot-wire anemometer, namely, (i) 
constant-current operation and (ii) constant- temperatu re operations. The 
constant current operation has the advantage of simpler electronic circuitry, 
however, constant temperature operation is more suited for measurements of 
fluctuations of high frequencies such as turbulent flows and measurements of a 
boundary layer where velocity varies from essentially zero near the wall to 
large value at the edge of the boundary layer in a small boundary layer thick- 
ness. Another advantage of constant temperature operation is the possibility 
of obtaining greater flow sensitivity using high overheating with no risk of 
probe burn out owing to sudden velocity decrease. For these reasons the constant 
temperature operation for hot wire anemometer was used during the present 
experimental measurements in the case of blunt trailing edge airfoil test. The 
constant temperature operation is briefly described in the following paragraph. 

Figure Ml-l6 shows the schematic diagram of the principle of constant 
temperature operation. The idea behind the constant temperature system is to 
minimize the effect of thermal inertia of the probe by keeping the sensitive 
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FIGURE HI-15 - RELATION BETWEEN ANEMOMETER VOLTAGE AND FLOW 

VELOCITY WITH AND WITHOUT THE USE OF L3BIAHIZBI 
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element at a constant- temperature and using the heating current as the measure 
of heat transfer and hence velocity. The principle requires a sophisticated and 
well-designed electronic system, and it can be explained with the aid of Figure 
111-16. Under conditions of bridge balance a voltage is present across the 
vertical bridge diagonal. This voltage is supplied by the servo amplifier. A 
slight change in the convective cooling of the sensor will cause a small voltage 
to appear across the horizontal diagonal. The latter voltage, after having 
undergone considerable amplification, is fed back to the vertical bridge dia- 
gonal, its polarity being selected so that it will automatically balance the 
bridge. In this way the temperature variations of the hot wire are kept 
extremely small. 

I I 1.3-3 Probe Dr ive Meehan i sm : Boundary layer surveys, both on the upper 

surface of the airfoil as well as in the wake behind airfoil trailing edge, were 
obtained by using an automatic probe drive unit to position the probe at pre- 
selected heights above the model surface. This unit, which is shown in Figure 
111-17, was mounted on the outside wall of the test section with the probe 
extending through the wall into the test section. Figure MI-18 shows the 
typical model installation in the wind tunnel facility and the boundary layer 
probe arrangement. This boundary layer probe can be located at any chordwise 
location on the airfoil as well as in the wake aft of the airfoil trailing edge. 
This probe is driven in the direction perpendicular to the local airfoil contour 
or normal to the airfoil chord line extended in the airfoil wake. The probe 
drive units consist of a 28-volt gear motor, lead screw, drive nut, and position 
potentiometer. A maximum of 50 pre-selected probe positions may be input to the 
data acquisition program. During the run the program will signal the gear motor 
to drive the probe assembly to each position in sequence and pause while the 
data are read by a scanivalve and transducer located outside the test section. 
Minimum increment of probe travel is .0318 centimeters. Maximum travel is 8.9 
centimeters. Limit switches are provided to protect the mechanism and prevent 
the probe from crashing into the model. 

A multi-post scanivalve and pressure transducer were used to read the probe 
pressures. The transducer is rated at 0.175 kilogram per square centimeter 
di fferent ial . 


|j|,4 Data Acquisition, Reduction and Analysis Technique 

Figure IM-19 shows the schematic diagram of the system used for the 
acquisition and reduction of the raw test data. The heart of the system is a 
Lockheed Electronics MAC 16 computer. The raw data was made available, in an 
abbreviated form, on teletype for outline monitoring of the test, and in its 
entirety on paper tape using a high speed punch. Final data reduction was 
accomplished on a Uni vac 1106 Central Computing System with remote access 
terminals. General description of Data Acquisition Unit and MAC 16 computer is 
given in Appendix C. In this section, however, main equations used in the 
reduction of data are given and the technique and the computer program for the 
analysis of data in the airfoil wake is described in brief. 
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I I 1.4.1 Data Reduction Equations ; Tunnel conditions acquired at each data 
point consisted of the dynamic pressure q m , the total pressure H m , and the total 
temperature T m . In addition the barometric pressure P a is obtained from a 
thumb-wheel switch input. These conditions are corrected to free stream condi- 
tion using standard procedure. 

An initial estimate of the free stream Mach number, 


'vn 

est “ 0.7 P a (1 + .000171 q m ) 

from which the compressible dynamic pressure is, 

q = q m (1 + K (2 - 

where K is the solid blockage factor calculated as, 


<est» 


(MI-6) 


(MI-7) 


7 c 


Kl x V w Yl + 0- i » M estV 0 - 5C Dm"\ 

+ \y^~ j\ -W-J 


(MI-8) 


where 


W w 

Cl 

Co 

m 

c 


and 


W 


* wing blockage shape factor 

t = blockage correction factor relating model span to tunnel width 
= model volume (cubic feet) 

= tunnel cross-sect ional area (square feet) 

= model drag coefficient 
= model chord 
= tunnel width (ft) 


The free stream total pressure is 


H m + \ % 2 * % 


(MI-9) 


where and K^ 2 are airstream calibration constants for the test section. 


The free stream static pressure, mach number, static temperature and 
density can now be calculated as, 
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T m + 459-fr 
00 1 + .2 M 2 ~ 


P 


00 


.0005826 



(I 1 1-10) 


The finally reduced data assumes Incompressible flow conditions and are there- 
fore based on ' ‘ 1 


^o = 


qco 


1 + 




+ 



1000 


(111-11) 


and 


Po - p. 0 - -2 M„ 2 ) 2,5 . 

Denoting the static pressure measured during the surveys as P 5 then, 


C P = 


P s • P 

5 a 


from which the local Mach number and density are obtained... 


M e = 


1 + 0.2 M 


(1 + .7 M^ 2 C p )- 286 


-1 


1/2 


(I I 1-12) 


(II 1-13) 


and 


- + H , 2 \ 2 - 5 

1 + .2 M e 2 / 

I I 1.4.2 Data Analysis Technique : A computer program was formulated for the 

data analysis of the experimental measurements obtained in the viscous flow in 
the wake. ■ Exper i menta 1 1 y measured variations across the airfoil wake, values of 
total pressure, static pressure, ambient temperature, free stream dynamic head 
and barometric pressures at two closely spaced X locations constitute the input 
to the computer program for the data analysis. This program then calculates 
and prints out the following information: 

(1) nond i mens i ona 1 velocity profiles at two stations, 

(2) shear stress distribution across the viscous layer between two 
stations , 
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(3) mixing length distribution, eddy viscosity ratio distribution, and 
slope of the velocity profiles at various y locations for the 
station between input stations X x and X 2 , 

(4) turbulent dissipation or shear work integral, displacement thickness, 
momentum thickness and energy thickness at both stations, 

(5) similarity parameters for velocity profile and values of similarity 
functions at both stations, and 

(6) various integrals at each station, such as 



where yi is any distance above the lower edge of wake in the viscous region and 
u( y ^ is the velocity at distance y above the lower edge of wake. 

In order to calculate the shear distribution, eddy viscosity ratio and 
variation of mixing length across the airfoil wake the use of the following 
basic equations was made: 


pv = " j ^ (p u ) d * + p o V o 
o 

Du Du Dp Pt 

Pu ax p v a Y a x * 


ay 


The Euler equation at the upper and lower edges of the wake is given by 


dP ei 

dx 


dlL 


dPe 


dUp 


- p e L U e|_ dx 5 dx " Pe u Ue u 


dx 


(II 1-14) 


(IN-15) 


(II 1-16) 


By integrating equation (l I 1-15) from y = 0 (lower edge of wake) to y = yi, 
making use of equations (111-14) and (111-16), and using Leibnitz's rule and 
integration by parts, the following equation can be derived after some 
algebraic simplification: 
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(MI-17) 


In the above equation the symbols have the following meaning: 

Yl = distance above the lower edge of wake 
U eL = average velocity at the lower edge of wake for stations X x and X 2 
u e u “ average velocity at the upper edge of wake for stations X and X 

“em " i <% * V 

y D = lower edge of wake 
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V Q = y component velocity at lower edge of wake 

P 6u = average value of pressure at the upper edge of wake 

P eL = average value of pressure at the lower edge of wake 

Pe u * average value of density on the upper edge of wake between 
stations Xj and X 2 

p e . = average value of density on the lower edge of wake between stations 
X x and X 2 

Pe m - i(Pe u + Pe L ) 

t( V i ) * value of shear stress at y = y 1 

t (YoV * va * ue shear stress at lower edge of wake 

* 0 

The right-hand side of equation (111-17) was programmed to compute shear 
stress T/y j at various distances, y^, above the lower edge of the wake. Values 
of shear' 1 stresses at the lower and upper edges of the wake were assumed to 
be zero. The'' same equation and the same numerical procedure was used to com- 
pute shear distribution, eddy viscosity ratio and turbulent dissipation in the 
various regions of the wake flow behind the airfoil trailing edge. 

I 1 1.4.3 Determination of Profile Drag from Experimental Measurements in 
the V/ake : The expression due to Betz, Reference 11, for the computation of 

profile drag of airfoil from the measurements of total pressures and static 
pressures across the wake behind the airfoil is given by: 



profile drag coefficient of an airfoil 

P T = total pressure of free stream 
'e 

static pressure across the wake 


where Cq = 
P T = 

00 

Ps (y) = 
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2 


q ro = free stream dynamic head = ipU^ 

P = free stream static pressure 

CO 

C = airfoil chord 

If static pressure, P s , * in the above equation (111-18) is assumed to be con- 
stant and equal to the'free stream static pressure, the equation (I 11-18) can 
be s impl if ied to, 

Co- (y/c , ' 

or 

C D 


or 


where 6^ = momentum thickness in the wake far downstream of airfoil trailing 
edge where static pressure has stabilized to the free stream value. 


= 2 




f | I .5 Types of Tests 

Experimental data were obtained for model under two conditions, namely (i) 
sharp trailing edge airfoil, and (ii) blunt trailing edge airfoil. In the case 
of blunt trailing edge airfoil, the trailing edge thickness was kept at 1 per- 
cent of the airfoil chord. Airfoil chord c is equal to 29.2 centimeters and 
the value of the maximum thickness of the airfoil tested is 15 percent. The 
chord of the blunt trailing edge airfoil is approximately equal to 0.975c where 
c is the chord of the airfoil with the sharp trailing edge. This is shown 
schematically in Figure 111-19 which also shows the coordinate system under 
which measurements are obtained during the tests. 

All the measurements taken during the tests were obtained at free stream 
Mach number of approximately 0.2 and free stream dynamic head of 60 psf. These 
■conditions would correspond to the free stream Reynolds number during the tests 
of approximately 1.3 to 1 . *» million based on the chord of the airfoil. Measure- 
ments were obtained at angles of attack between 0° and angle of attack corres- 
ponding to the occurrence of incipient separation on the upper surface of the 
airfoil in the neighborhood of the trailing edge; the latter condition corres- 
ponds to a - 10.8° for the test airfoil. Measurements taken during the tests 
can be classified by two categories namely (i) measurements of pressure dis- 
tribution on the surface of the airfoil and (ii) viscous flow measurements. 
Measurements of the viscous flow were obtained both on the upper surface of the 
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airfoil in the neighborhood of the trailing edge as well as in the airfoil wake 
at several chordwise locations. 

In the case of test airfoil with the sharp trailing edge the viscous flow 
measurements were obtained by the use of total-static pressure probe which is 
shown in Figure MI-9. For the test airfoil with the blunt trailing edge, the 
viscous flow measurements were obtained by use of hot wire anemometer probe 
which is shown in Figure 111-14. Figures 111-20, 111-21 and 111-22 show the 
plots of measured pressure distributions and loci of minimum velocity for the 
airfoil with the sharp trailing edge. Boundary layer surveys in the airfoil 
wake have been obtained at x'/c locations shown in Figures 111-20, 111-21 and 
111-22 by triangular symbols. These measurements are tabulated in Table 3~3 
for sharp trailing edge airfoil. As seen in this table these measurements at a 
particular chordwise location x'/c in the wake consist of values of (i) total 
pressure P^, (ii) static pressure P s , (iii) free stream dynamic head q^, and 
(iv) ratio of u/U^ at several y locations in the wake flow boundary layer. 

Figure 111-23, 111-24 and 111-25 show the plot of measured pressure distri- 
butions on the airfoil surface and the loci of minimum velocity in the wake of 
blunt trailing edge test airfoil. Boundary layer surveys are obtained, by the 
use of hot wire anemometer probe, at the chordwise locations s'/c indicated by 
the triangular symbols on measured loci in the above figures. The boundary 
layer surveys consist of measurements of u/U ro or u/U e and q TO for several y 
values and for a given chordwise location x'/c. These measurements are 
tabulated in Table 3~4. 


IV. RESULTS AND DISCUSSION 


In the previous section the description of experimental facility, type and 
number of experimental measurements and the computer program for the analysis of 
the fundamental wake flow parameters was given. From the output of this computer 
program various physical parameters for the flow in the wake of single component 
airfoils are studied. The relationships among various physical parameters, which 
appeared in theoretical equations of Section II, are derived and described in 
this section. Establishment of these relationships between various physical 
parameters is of vital importance for the prediction of drag of arbitrary con- 
figuration single or multi -component airfoils by the generalized theoretical 
method such as described in this report. 


IV. 1 Presentation of Typical Measured Experimental Data 

Figures I V— 1 through K V~5 show the measurements of velocity profile on the 
upper surface of the test airfoil with sharp trailing edge. These measurements 
were performed by use of total-static pressure probe combination shown sche- 
matically in Figure I I I ~9 • This type of boundary layer measurements on the 
surface of the airfoil were obtained for the purpose of determining (i) the 
approximate location of the point of incipient separation on the airfoil 
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Table 3~3 - MEASUREMENTS IN THE WAKE OF SHARP T.E. TEST AIRFOIL 
THE USE OF TOTAL STATIC PRESSURE PROBE COMBINATION 
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Table 3~3 (Continued) 
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-.05** 

*1.54 

*112.5 

7052-9 

.90*4 

-.o«o 

41.90 

*111.5 

*052.1 

.9094 

-.0449 

*1.4* 

2111.0 

2052.9 

.97<* 

-.0495 

<1.94 

*11*. 35 

'205*. 1 

.9001 

-.0*17 

<1.1* 

2111.0 

*055-5 

•9**7 

-.05*2 

41.24 

2111.3 

2052.1 

.9*11 

-0*5 

*1.90 

*10*.4 

2055-0 

.944* 

-0*<9 

41.47 

*110.0 

0052.3 

.9470 

-.OJ06 

<1.10 

2107.2 

*053-05 

.9294 

-0117 

41.4* 

*100.75 

*052.3 

.9547 

-.0*0 

<1.*7 

2104.2 

*054.4 

.0979 

-.0345 

<1 .23 

*107.3 

*05*.* 

.9414 

-.0*17 

<1.21 

*101.1 

>054.< 

.•<*5 

-.090* 

41.35 

2105.25 

0055.25 

•91T7 

-.0174 

<1.15 

*040.4 

2054.7 

.•410 

-.0040 

(1.1* 

*103.0 

1053.5 

.90H 

-.01 » 

<1.4* 

2-9<.4 

2055.0 

.*109 

-.0*17 

41.29 

2101.35 

*053.75 

.0041 

-.00*7 

<1.44 

*>94.7 

2054.9 

.0004 

-.0174 

41.00 

2101.2 

■054.0 

.474* 

-.00904 

<1.* 

*094.7 

2054.7 

.4044 

-.0139 

<1.40 

*100.5 

*055.95 

.out 

.00*1 

<1.15 

20*.4 

2054.95 

-4*51 

-0047 

<1.54 

*099.45 

*053.9 

.■<07 

.007* 

<1.*7 

*09*.7 

*054.1 

.•49* 

-0304 

41.49 

*099.4 

*053.0 

.0593 

.01*1 

<1.*9 

aiot.5 

*053.7 

.0490 

.00441 

(1.55 

*100.* 

2053.4 

.•4*7 

.0174 

<1.* 

*105.7 

0053-4 

.9205 

-0074* 

<1.4* 

*H».*5 

*095.9 

•»?94 

.0*17 

<1.50 

*10?. • 

0053.* 

•9404 

.01*1 

<1.51 

2102.25 

*055.*5 

.•909 

.0*1 

<1.55 

*104.4 

*05*. 9 

•95<7 

.0174 

41.47 

*KH.O 

*055.0 

-90*» 

.0504 

<1.49 

*110.4 

*052. • 

.94*0 

.0*17 

41.4* 

2105.4 

0052.5 

.9274 

.0540 

<1.49 

2111.4 

*05*.* 

->777 

.0041 

41.3* 

*107.15 

0052.25 

•9441 

.099* 

<1.44 

*112.1 

*051.3 

.9*27 

.0304 

41.54 

2100.4 

2052.1 

.9549 

.0495 

<1.* 

2112.2 

2052.5 

• 90f5 

.0944 

41.51 

*W9.5 

*052.0 

.9455 

.0*1 

<1.55 

*11*. 5 

2052.5 

• 9*43 

•091* 

41.40 

2110.* 

*051.0 

.9712 

.0*09 

<1.20 

2112.5 

2052.4 

.9*4* 

•0435 

41.55 

2110.05 

*051.9 

.9775 

.0*95 

<1.* 

2112.5 

2051.4 

■9*40 

.05*1 

41.47 

2111.4 

1051.# 

♦90*9 

.0*95 

<1.* 

2112.5 

*052.4 

.9*4* 

*0<09 

41.41 

*111.9 

*051.7 

•9*7* 






.04*5 

41.50 

2112.5 

0091.4 

.9954 






.074* 

41.40 

2112.3 

0051.4 

.9954 


Station 7» kl!o«rW»<juar* Mt«r 

*Vc ■ 1.107| a - 0.0* 

Praa Straw Total Praaaura/o f - 211 2. 5 
Praa Straw StaUc Praa«ira/C_ f - 20*7.9 
rni Straw Total Taaparatura ■ 305* * 


y/o 

Saf 


<raf 

. c 

-.0905 

(1.45 

2112.5 

2052.7 

•9**4 

-.0*00 

(1.33 

2111.5 

2052.7 

•mi 

-.0495 

(1.37 

2112.3 

2052.9 

.9*51 

-.050* 

41 .44 

2111.4 

2052.95 

.9772 

-.0444 

41.31 

2110.2 

2053.1 

.965* 

-.0417 

41.57 

2109.1 

2053.5 

.95*9 

-.0345 

<1.35 

*107.* 

2053.7 

•9399 

-.0304 

41.3* 

2106.5 

2055.* 

.9*77 

-.0240 

4i.3« 

2105.9 

2055.9 

.9141 

-.0*17 

41.54 

2104.05 

2054-0 

.9040 

-.0174 

41.29 

2103.2 

*054.115 

•*949 

-0t>9 

41.45 

2102,7 

2054.2 

•**9* 

-.0067 

61.2* 

2102.3 

*054.2 

.M62 

-.00904 

41.21 

2101.7 

*054.1 

.*•14 

.002(1 

61.41 

2101.6 

2055-9 

•***5 

.007** 

61.40 

2102.5 

2055.7 

.*927 

.0121 

41.44 

2053.6 

2055-6 

.*9*2 

.0174 

61.35 

2103.9 

2053.55 

-906a 

.0217 

61.41 

2104-5 

2053.4 

•9136 

.02(1 

61.52 

2105.5 

2055.* 

-9243 

-0304 

61.49 

2106.7 

*055-1 

•9357 

.034* 

61.42 

2107.7 

2052.9 

• 946* 

.0592 

61.55 

210*. 9 

2052.7 

• 95*2 

-0435 

61.32 

2109.45 

2052.6 

.963* 

.0521 

61.32 

2110.7 

2052.5 

.975* 

.0609 

61.3* 

2111.7* 

2052.3 

• 9*59 

.0693 

61.19 

2112.5 

2052.3 

•9916 

.07*2 

61.26 

2112.5 

2052.3 

.991* 


Table 3*3 (Continued) 


*8 

,&V 

* 4.88 Kilofran/aquart Malar 
a. 8.o* 

Station Zf q 
*’/e - 0.223?* 

• *.86 Kllofran/aquara wtar 
• 3 8.0* 

Station 3 j 0 , 

*v* - o.jur' 

• 4.88 nuingFapaia Malar 
• - 8.0* 

Fr*a Strata Total Fraaawra/o f - 
|»ra* Straaa Static Fraaaura/q" f - 
Fraa Straaa Total TaMoaratura « 

2120.0 

2058.1 

Fraa Straw Total Praaaura/q . - 
Fraa strum Static Fraaaura/S",, « 

2120.0 

2058.1 

Fraa Straw* Total Fraaaara/q^, - fctvf.O 
Fra* Straw StaUa fiawM/a^P*. • 80*0.8 

* 

Fraa Stria* Wtal Taaparaiura” 1 • 

306.1* * 

fra* straw Total Taaparatwra^ - 304.8* ■ 

Fra* Tranaition 



Fraa TranalUt* 




Fraa TrwaiUan 




F/t 

S*t 

41.3* 

r t 

Vaf 

> t 

«raf 


7/c 

Sat 

r t 

Saf 

f 

«nf 

* 

t/« 


r. 

A. 

w 

<T 

.30*0 

*1*0.0 

*047.7 

.9340 

.2345 

41.34 

*117,0 

2041.5 

.5484 

.8488 

414» 

*117.0 

*848»3 

*))■• 

.3175 

41.30 

8118.1 

*048. t 

.9890 

.25*0 

i.<U54 

*117.0 

1041.5 

.9484 

.8498 

41.18 

8117.0 

8840.3 

.9349 

.387* 

<1.07 

8118.0 

8049.2 

.8390 

.*<93 

<1.28 

*117.0 

*041.5 

.9481 

•8*18 

41^1 

*117.8 

8040.9 


>3380 

<1.88 

8998.4 

*070.1 

.4480 

.2*70 

41.79 

*117.0 

*041.3 

•94*4 

.3048 

41.fl 

*117.8 

8*80.9 

•9I88 

3)7« 

<1.11 

*084.9 

*070.1 

.5*50 

.304* 

<1.*3 

*114.9 

8041.5 

.9484 

•18*8 

41.18 

8114.4 

8808.19 

•M 

344* 

<1.15 

8IMB.1 

*-49.9 

■ .*400 

.3130 

41.34 

2114.7 

2041. r 

.9480 

•390* 

(41.19 

8115.8 


33«> 

<0.9< 

*089.2 

*048.8 

',4880 .V 

.5*80 

41.17 

2114.1 

80(1.0 

■9459 

•3390 

41.81 

*113.3 

8840.1 

•00 

3540 

<1.11 

*087,9 

*047.5 

.5790 

.3307 

41.27 

*114.85 

8041.0 

■9MJ 

•34*0 

41.91 

81*9.8 

8848^ 

•88K 

.3410 

<1.10 

8094.3 

*M.i 

.4800 

.3390 

4050 

8W9.3 

8040.5 

.8*13 

•3370 

<1.31 

8101.8 

8848.0 

.811* 

.3440 

41.14 

*100.5 

8045.? 

.7*80 

.3440 

41.54 

*103.73 

8040.4 

••984 


41.43 

8098.9 

88)9.9 

.791* 

.37*0 

41.09 

1108.0 

*045.7 

.8340 

.34*0 

<1.35 

8057.5 

8040.83 

•777* 

.3490 

41.84 

•*93.7 

•099.81 

.740 

.3773 

<1.13 

*115.8 

2044.9 


.35*0 

41.25 

7098.1 

80(0.0 

>7*84 

•3<95 

41.89 

8094.* 

•899.9 

.9444 

.3*84 

<1.81 

*117.2 

2047.J 

.9*40 

.3547 

*1.15 

8089.1 

80(0.0 

.4881 

•JT98 

41.18 

8099.5 

•*98.8 

*im 

.3*00 

<1.05 

*119.1 

*047.8 

.9*10 

.9400 

41.33 

8088.0 

8040.2 

.47*4 

•jm 

41.30 

1093.4 

8099,8 

.•884 

.9340 

41.00 

2119.8 

2048.4 

.9*oo 

.3450 

41.45 

8085.5 

8040.3 

.(881 

.3*10 

41.19 

ae.1 

88)1.* 

.71*4 

.4075 

<1.15 

21*0.0 

*048.4 

.9*80 

.3<95 

41.41 

8058.8 

80(0.5 

.7*37 

.3*70 

41.88 

80VT.8 

8099.73 





.3740 

41.45 

*0*4.5 

MO. 4 

• 7489 

.39» 

<1.75 

tin.) 

8898.7 

.00 






.5*0 

41,35 

*ioj.3 

20(0.8 

*8380 

.4000 

41.89 

MM.) 

8*98.4 

.804 






*3900 

41.45 

*110.0 

80(1.0 

.891* 


41-33 

8118,4 

8898.9 

.9*8* 






.3995 

<1.*9 

*113.8 

20(1.5 

-9398 

.*110 

41.18 

*114.4 

8898-4 

.940 






.4080 

41.30 

21K.4 

80(1.5 

.9501 

•4840 

41.89 

8115.8 

W98-4 

.00 






•4«3 

41.19 

*t!7.0 

*0(1.2 

.9510 

.4998 

<1.81 

81*4.3 

80)8.9 

.00 






♦4*40 

41.51 

2117.0 

20*1.0 

* .9587 

.4431 

(0.88 

*117.0 

8*98.83 

.00 






.4340 

<1.50 

*117.0 

80(1.0 

.95*7 

.4518 

41.1* 

*117.0 

*859-*9 „ 

.990 






.4430 

<1.21 

*117.0 

80(1.0 

.9587 

•4410 

41.88 

8117.0 

8859.13 

.00 


Station * i q . - *.90 tlloarw/aquara aatar 
x'/c - 0.517™ 

Praa Straw Total Praaaura/q . « 21 *7.° 
Fra* Straw Static Praaaura/6*~ f » 2050.8 
Prat Straw TVjtal Taaparatura • J08* K 
Fraa Tranaition 


r/o 

5sl- 

W 

*t 

P 

«r.f 

* 

•*T55 

(1.40 

*117.0 

20(0.3 

.9581 

.29*5 

41-34 

2D7.0 

8060.5 

.95*4 

.3100 

<1.35 

2117.0 

2060.3 

.950 

.5*75 

<1.27 

2114.7 

2059.9 

.9578 

.5450 

41.53 

2114-4 

2059.6 

.941* 

.5540 

(1.35 

2il*.3 

2059.3 

.9*3* 

.3*25 

*1.40 

2108.5 

20W.5 

.*#» 

.37»0 

*1.25 

*108.2 

2059.45 

.0313 

.3800 

*1.51 

*099-4 

2039.4 

.•041 

.3840 

61.35 

3098.8 

2059.55 

.7984 

.3*95 

61.30 

2099.0 

2059.3 

.8011 

.3930 

61.04 

2099-5 

2059.25 

.80(7 

.3973 

(1.70 

*100.0 

2059.2 

.8112 

.4015 

61.25 

*101.0 

2059.2 

.8*21 

.40(9 

61.26 

2102.X 

2059-15 

• 8343 

.4100 

61.22 

2104.0 

2059.15 

.8513 

.4130 

61.28 

2104.7 

2059.15 

.8768 

.4242 

61.26 

2111.0 

2059.1 

.31(0 

*4380 

61.17 

2113.85 

2059.0 * 

.5417 

.4415 

61.33 

2115.3 

2059.0 

.9558 

.4500 

61.23 

2116.5 

2059.0 

.9651 

.4590 

61.14 

2116.9 

2059*0 

■9693 

-44TO 

61.17 

2117.O 

2058.8 

•9705 

.4750 

61. 31 

2117.0 

2058*75 

-9705 

.4*50 

61.21 

2117.0 

2068.75 

.9705 


Station 5j q„, - *.88 813.o#rw /apiara «atar 

«'/» ■ 0.885?* f a ■ 8.0* 

Fro* stna Total Fraanira/d • 2117.0 
Fraa straw Static FmaauraA" - *05*.* 
FW Straw Total T*a*»r*t»r* - 312.«* * 

Fro* Tranaition 


fA 


* 

r 

f 

Vaf 

S,f 

V.f 

w aa 

.5440 

61.10 

8117*0 

80*8.03 

■ 9*10 

• 9570 

61.84 

1117*0 

*>5*.03 

.9*W 

.3650 

61.29 

*117.0 

*030.01 

.9810 

• 5740 

61.15 

2117.0 

8098.* 

.9808 

.909O 

<1.81 

*115.5 

8058-1 

.940 

.5915 

(1.13 

H15.5 

8050.15 

.95» 

.4000 

(1.18 

8110.9 

80)8.19 

•9098 

,4086 

<1.33 

*108.3 

>058.13 

.900 

.4175 

<1.14 

*106.3 

2058.1 

.8881 

.4*60 

<1.90 

2*09,0 

2058.9 

.8700 

.4350 

<1.0 

2M4.3 

K69.5 

.8)19 

.4490 

(l.lf 

2104.3 

80)9.7 

.8)40 

.4515 

61. 14 

*104.8 

1099.7 

.8)80 

.4<10 

<1.10 

8)84.1 

*80.5 

.*781 

.4700 

(1.19 

1107.9 

800.3 

•001 

•4785 

41.10 

1109.8 

*0)9.3 

.9041 

«48|0 

41.21 

2111.8 

*059.0 

.980 

.4950 

(t-07 

1H3,4 

*0)0.5 

.9443 

.5052 

61.0J 

8114.8 

*050-3 

.960 

.5190 

61.09 

81 19.4 

1050.1 

.9(08 

.52*0 

61.20 

8117.0 

8067.9 

-900 

.5900 

61.81 

8117.0 

80)7.9 

.9*90 

.3390 

61.81 

81 17.0 

*057.8 

.98)0 



Table 3"3 (Continued) 


Statlw ij a . - *.88 ri 1 n|r w'ty nur 

*'/e - .0**ta- io.8- 

Praa Straw Total Praanua/a - - 2«2i.) 

Praa Straw Statia Pracaura/l”- . 2060.1 
Praa Straw Total T^mrmturm 1 - 311.1" X 
Praa Trwaitiw 


luua 2, . . 6.88 uim^tptn attar 

«V© - .ofclT e- 10.®* 

Tital PhibiVl , ml 
IU 1 K - 1 

Total tapr«t«i Nr • J 


Praa Itna Total Praanua/t 
Tr— ItM “ “ 

8tm Stna 


. 2111 .) 

'1.15 
» 5»i.»* x 


x/o 

Saf 

Sar 

P 

Saf 

Ht 

*/• 

-2*. 

Sar 

P t 

S*f 

Saf 

* 

..Otfi 

61.45 

2121.5 

206*.3 

.9)04 

-.0061 

61.4) 

*1*1.) 

*060.) 

•90* 

-.0*17 

61.1) 

21*1.5 

2068.3 

.5304 

-.0*17 

61*08 

21*1.) 

*868.) 

•9*6 

-.0174 

61.85 

.21*1.5 

S068.3 

•9)04 

-.0174 

61.46 

*119.0 

*060.6 

.9098 

-.01)1 

61.37 

*1*1.0 

*068.4 

.9*31 

-.01)0 

61.46 

•117.6 

*068.6 

.4*8 

-.01)1 

61.45 

* 119.7 

*06*. 5 

.91*7 

-.01)1 

*1.1* 

tilt.) 

*06*-6 

•831 

-.0109 

61.4* 

* 117.4 

*068.6 

.*» 

-.OWf 

61.47 

nou 

8068.7 

.** 

-.00*7 

61.36 

X11X.5 

*068.4 

.0451 

-.00*7 

61.81 

* 100 .* 

*0*0.7 

.7*0 

-.0065 

61.3* 

2104.3 

*068.7 

.76)1 

-.006) 

61.4) 

•094-* 

*NM 

.<4* 

-.0045 

61.4* 

*056.0 

*068.4 

.64)5 

-.004) 

61. St 

*0*8.7 

809.* 

■9691 

-.008* 

61.45 

*0*7.5 

*064.* 

•3455 

-.00*0 

61-40 

*0*4.4 

**9-4 

.497* 

0.0 

61.35 

20*0.* 

»65.5 

.42*7 

0.0 

61. )T 

*0*1.) 

*69-3 

.48* 

.008* 

61.40 

2076.0 

*06).9 

.3145 

.00*0 

61.41 

*019.0 

*69.) 

.89* 

.0044 

61.44 

*074.1 

*010.2 

.2)10 

.0044 

61.JJ 

*077-0 

MM 

.869 

.0065 

61.41 

*073.* 

»70.5 

.*055 

.006) 

61.44 

00TT.5 

8069.8 

.89* 

.00*7 

61.33 

*075.0 

*070.4 

.1056 

.00*7 

61.J) 

*07*.© 

0090.0 

•JM 

.0131 

61.56 

*074.3 

*070.1 

.*61) 

.01)1 

61.J9 

*0|*.7 

8070.0 

•sm 

.0174 

61.47 

*076.3 

*065.9 

.3*76 

.0174 

41.J1 

*0*8.5 

*069.4 

-4544 

.0817 

61.41 

*075.4 

*065.5 

.401* 

.0*17 

*1*81 

*0*5-7 

*69.0 

.90* 

.0861 

61.55 

*0*5.0 

*069.4 

.470) 

.0*61 

41.88 

*0*9.) 

*069.7 

• 9694 

.0)04 

61.54 

*0*7.1 

*064.) 

.5)59 

.0)04 

*1-87 

*096.0 

*69.3 

.6)1« 

.0)47 

61.80 

*0)1.7 

3069.* 

.6049 

.0)47 

41.)* 

*099.* 

*069.5 

.49* 

.09)1 

61. 51 

*047-* 

*064.1 

.6760 


41.47 

1M4-7 

*69.4 

.19* 

.0434 

41.45 

2108.* 

*065.1 

.740) 

.04)4 

41.)* 

2105.* 

*69.0 

.4*6 

.0477 

61.50 

*108.0 

2064.0 

.7564 

.0477 

41.44 

*109.) 

*60.7 

.79* 

.0)** 

61.54 

*111.3 

*064.0 

.*411 

.0)80 

41.4) 

*W)-9 

*06*.) 

.70* 

.0) 6) 

41.54 

*115.) 

*068.8 
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41.2) 
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41.24 
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41.24 
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.9750 

*1.95 
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-.00)8 
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41.80 
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.*950 
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• 92)0 
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41.2) 
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41.14 
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41.92 
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.0174 
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41.05 
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41.1) 
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41.14 
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.0*95 
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41.1* 

•0 *95 
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*1.09 
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.9*4* 
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.9*40 
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41.23 

.1477 
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.0414 .9)00 61.5* 
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.045) .9490 41.29 
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41.1* 
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.7740 
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.0015 
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-.0331 
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-OW* 
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.71® 

*1.30 
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-.003* 
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*1.37 
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*1.33 

.0113 
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41.24 
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.4124 

.03** 
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FIGURE III-20 - MEASURED PRESSURE DISTRIBUTION ON AIRFOIL SURFACE 

AND LOCUS OF MINIMUM VELOCITY AT 0< - 0° 

(FREE TRANSITION) 
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FIGURE III-22 - MEASURED PRESSURE DISTRIBUTION ON AIRFOIL SURFACE 
AND LOCUS OF MINIMUM VELOCITY AT oC- 10.8* 
(FREE TRANSITION) 
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FIGURE III-2J - MEASURED PRESSURE DISTRIBUTION ON AIRFOIL SURFACE 

AND LOCI OF MINIMUM VELOCITY AT c< - 0° 

FOR BLUNT TRAILING EDGE AIRFOIL (FREE TRANSITION) 




BLUNT TRAILING EDGE AIRFOIL AT C<« 4.0* 
(FREE TRANSITION) 
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FIGURE III-25 - MEASURED PRESSURE DISTRIBUTION ON AIRFOIL SURFACE 
AND LOCUS OF MINIMUM TELOCITY AT cX. * 10.7* for 
BLUNT TRAILING EDGE AIRFOIL (FREE TRANSITION) 




surface and, ( i I ) to obtain the initial ccmd it ions for the purpose of wake flow 
calculations by the theoretical profile drag method developed under the present 
studies. Figures I V— 1 through IV^4 show the results of measurements of 
velocity profiles on the upper surface of the airfoil for angles of attack a of 
0°, 8°, 11° and 13°, respectively ; ohqrdwise Vocations of measurements are also 
indicated in these figures. It can be seen frojn these figures that, for angles 
of attack of 0° and 8° and for X/C locations 'forward of or equal to 0.9 in case 
of a of 8°, the shape of velocity profile is smooth and there is an absence of 
scatter in the measured values.^ However, for angles of attack a of 11° and 13°> 
and in particular for chordwise location close to the trailing edge, the shape 
of the velocity profile is not smooth near the surface as evidenced by the 
appreciable scatter in the measured data. The cause of the above phenomena is 
found in the ability of total-static pressure probe combination . to measure 
reliable velocity profiles under various flow conditions. The pressure probe 
of Figure I II —9 is reliable when the viscous flow measurements are made for 
attached flow conditions. However, the measurements of velocity profile by 
this probe, for the boundary layer which is approaching separation or in 
separated flow region, are not reliable and are subject to interpretations. 

Hence, the measurements by this probe design under the above-mentioned conditions 
can be used for the purpose of qual i tat i ve analys i s only. However, in order to 
determine the position of initial separation and the extent of separated flow 
region, measurements performed by the present design of the pressure probe can be 
effectively used for this purpose by the procedure such as shown in Figure I V— 5 - 
Various investigators have observed from experimental measurements of velocity 
profiles that the value of the form factor H at initial separation point for 
ordinary turbulent boundary layer is greater than approximately 2.0 and that 
exact value of H at separation depends upon surface roughness. Thus, in order to 
determine the initial separation as well as the extent of the separation region 
from pressure probe measurements, growth rate of integral thickness, such as the 
displacement or momentum thickness, needs to be examined in addition to Form 
Factor H. This is done in Figure IV-5 and the following conclusion can be drawn 
from this figure: 

(1) Separation at T.E. starts at ot = 8.0° 

(2) Separation at X/C = 0.9 starts at a = 11.5 to 12° 

(3) Separation at X/C = 0.8 starts at a = 13*0° 

Figures IV-6, IV- 7 and IV-8 show the plots of velocity profile measured by 

the hot wire anemometer technique. These velocity profiles were measured on the 

upper surface of the airfoil with the blunt trailing edge of 1 percent thicknesses. 
Angles of attac.% corresponding to Figures IV-6, IV-7 and IV-8 are 4°, 10.8° and 13° 
respect i vel, . It can be seen from these figures that measurements of velocity 
profiles by hot wire anemometer are smooth and reliable for attached boundary 
layer flow as well as for viscous flpw approaching separation. Hot wire anemometer 
technique is unable to give any i nd ication abput the direction of the flow in the 
boundary layer; i.e., it gives indication only of the absolute magnitude of 
velocity. In addition, as the magnitude of the flow velocity approaches zero at a 
point in the flow, the relative magnitude of the velocity of fluctuation is large 
compared to the mean flow velocity at that point. Under this circumstance the 
accuracy of measurements by the hot wire anemometer techniques suffer as low 
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FIGURE IV-1 - VELOCITY PROFILE MEASUREMENTS ON UPPER SURFACE OF 
AIRFOIL AT c< - 0* AND FOR SHARP TRAILING EDGE 
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values of velocities are approached in a turbulent flow with reverse flow veloci- 
ty profiles. Because of the above two reasons, it is necessary to exercise some 
judgment in fairing the velocity profile shape when measurements are made with 
the hot wire anemometer techniques in a boundary layer flow having reverse flow 
velocity profiles. Moreover, large numbers, of data points, in the vicinity of 
values of y locations where the zero velocity is approached, would make it 
possible to exercise sounder judgment in fairing reverse flow velocity profiles 
near the y locations when the values of velocity are near zero and near the wall 
where negative velocities are encountered. An example of this type of fairing 
is shown in Figure IV-8 for value of X/C = 0.965- 

An example of the type variations of velocity profile in the wake of the 
test airfoil with sharp trailing edge is shown in Figure IV-9- These wake 
surveys correspond to an angle of attack a of 8° and when transition on the 
lower and upper surfaces of the airfoil were fixed at 10 percent chord by the 
use of grit strips. The wake surveys shown in this figure were performed by the 
use of total-static pressure probe combination shown in Figure III-9- Near the 
trailing edge the rate of increase in minimum velocity with respect to the 
distance X'/C is quite large compared to the large distance downstream from the 
trailing edge. This figure also points out that the rate of variations in the 
width of the wake is small at a far distance downstream compared to near the 
trailing edge. It should be mentioned that this figure is not plotted to give 
true relative y locations for velocity profiles at different X*/C locations but 
to indicate relative shape of the velocity profiles at different X'/C locations 
in the wake behind the airfoil. 

Figure IV-10 shows the example of experimentally measured variations of 
momentum thickness and form factor in the wake behind the sharp trailing edge 
test airfoil. The data shown in this figure are for an angle of attack of 8°, 
for conditions of free transition on the airfoil surfaces and for free stream 
Reynolds number of approximately 1.3 million. It can be seen from this figure 
that values of momentum thickness and form factor decrease in the downstream 
direction in the airfoil wake and at far distance downstream in the wake their 
values become approximately constant. 

Examples of typical measurements of static pressure distributions across 
the airfoil wake behind the sharp trailing edge test airfoil are shown in 
Figures I V— 11, IV-12 and I V- 13; corresponding measurements of velocity profiles 
are also shown in these figures. Figure IV-11 is for angle of attack a = 0.0°, 
IV-12 is for a = 8.0° and fixed transition on upper and lower surfaces of the 
airfoil; and Figure IV-13 is for a = 10.8°. Measurements shown in these 
figures indicate that static pressure is not constant across the wake but has 
an approximate parabolic variation. These figures indicate that static pressure 
has a higher value at lower edge of wake and lower value at the upper edge of 
wake for angles of attack greater than zero degree. In addition, it is seen 
from these figures that the maximum value of the static pressure in the wake 
occurs in the vicinity of y location corresponding to minimum velocity point in 
the wake velocity profile. 

Figures IV-14, I V- 1 5 and I V - 1 6 show indirect measurements of shear stress 
profiles at angles of attack of 0° and 10.8° and at few chordwise locations in 
the wake behind the airfoil trailing edge. These shear stress profiles are 
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FIGURE IV-16 - COMPUTED SHEAR STRESS DISTRIBUTION AND EXPERIMENTAL VELOCITY PROFILE 
IN AIRFOIL WAKE AT ©< - 10.8° (X^C - 0.344. Xg/C - 0.517) 


computed by the numerical solution of integro-differential equation (111-17) in 
which the boundary conditions of experimentally measured velocity profiles and 
static pressure distribution is used. These figures show some interesting 
results such as the value of the shear stress T(y) is maximum (or minimum) at 
the y location in the wake velocity profile where the velocity has minimum 
value. Classical theories such as Boussinesq's eddy viscosity concept and 
Prandtl's mixing length theory would, on the other hand, indicate a value of 
zero shear stress on the locus of minimum velocity point. As a matter of fact, 
the value of shear stress, at a y location corresponding to the locus of mini- 
mum velocity, can become as high as 8 to 10 percent of free stream dynamic 
head as can be seen from these figures. The value of the shear stress on the 
minimum velocity locus is used as one of the boundary conditions in the 
theoretical equations in the present method for the calculations of the wake 
flow and the profile drag of a given airfoil. Hence, the use of Boussinesq's 
or Prandtl's hypothesis in the present method may lead to erroneous calculation 
resul ts. 


IV. 2 Investigations of Fundamental Parameters for Wake Flow 


{Theoretical equations in Section M contain such terms as (i) 


f(n) dn. 


f 2 (n) dn, 


(ii)' shear integral terms /x/p dy and shear distribution x(y 8 )/p on the locus of 
minimum velocity and (iii) integral for pressure distribution across the wake 
JTp dy, distribution of static pressure on the loci for outer and inner edges of 
wake and also on the locus of the minimum velocity. In order to solve these 
theoretical equations for the calculations of flow in the wake and prediction of 
profile drag of airfoil sections, it is necessary to express the above items by 
auxiliary equations as functions of either dependent variables or in terms of 
initial conditions in the vicinity of trailing edge. In this way the number of 
dependent variables, which are required for the complete specification of flow 
in the wake, are kept the same as the number of theoretical equations. These 
auxiliary equations are derived, as shown in the following paragraphs, by the 
use of dimensional analysis, investigation of generalized parameters and by the 
use of available results for "similar flows 11 . For this purpose use is also made 
of experimental measurements obtained during the present investigation. 


IV. 2.1 Functional Relationships for Velocity Profile Similarity for Wake 
Flow: Figures IV- 1 7 and I V- l"S are the plots of experimental data for similarity 

of” velocity profiles with different similarity parameters. Experimental data 
used in constructing these similarity curves are for chordwise locations in the 
wake which are far from the airfoil trailing edge. The same experimental data 
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FIGURE 17-17 - VELOCITY PROFILE SIMILARITY (l) IN REGION IV 
(FAR FROM AIRFOIL TRAILING EDGE) 
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for velocity profiles are used in constructing similarity curves of Figures IV-17 
and I V- 1 8 . Comparison of velocity profile similarity plots of Figures IV-17 and 
IV— 1 8 indicate that, when similarity parameter shown in Figure 1 V- 1 8 is used for 
the purpose of obtaining velocity profile similarity there is less scatter in 
experimental data than by the use of similarity parameter of Figure IV-17- This 
fact suggests that similarity parameter of Figure IV-17 is more suitable for true 
representation of wake flow characteristics at large distance downstream in the 
wake of single component airfoil with sharp trailing edge. For the purpose of 
making use of this information about velocity profile similarity, it is necessary 
to obtain an analytical expression for the relation between similarity parameter 
n and similarity function f(q)- The curve which fits the experimental data of 
Figure IV-17 has the analytical expression of the form given by 

1 0 884 

f( n ) = J [' + Cos (n a * * ir)]. OV-1) 

An analytical expression for the curve, which fits the experimental data of 
Figure IV— 1 8 , is given as 

1 0* 885 - 

f (n 2 ) 2 t 1 + Cos (n 2 ■ j)]- ( l v— 2 ) 

Equation ( I V— 1 ) and (IV-2) show the functional relationships for the 
similarity parameters for the viscous flow in the wake at chordwise locations far 
away from the airfoil trailing edge. Similarity plots for velocity profiles for 
the viscous flow in the wake of airfoil with sharp trailing edge and in the 
vicinity of airfoil trailing edge are shown in Figures I V- 1 9 and IV-20. Defini- 
tions for similarity parameters and variables in Figures I V- 1 9 and IV-20 are 
similar to those of Figures IV-17 and IV-18, respectively. As seen in Figure 
IV-19, the scatter of experimental data for wake flow in the neighborhood of 
trailing edge is larger than corresponding plot for far distance wake flow shown 
in Figure IV-17; the similar remark can also be made for data shown in Figure 
IV-18 and IV-20 which are for large and small distances in the wake behind the 
airfoil trailing edge. An analytical expression for the curve fit of the mean 
curve shown in Figure IV-19 is given by 

i 0- 825 ~ 

^ (n 3 ) = J [1 + Cos (n 3 • ^)] (IV-3) 

As the scatter in experimental data, shown in Figure IV-20 is larger than 
corresponding similarity plots of Figure IV-18, three different curve fits are 
defined by an expression of the following form: 






f (n, 

.) * 

J [1.0 + Cos (n 4 

where , 

for 

the 

curve 

Cl, 

p = 

1 . 0 , 


for 

the 

curve 

c 2 f 

p = 

0.755, 

and 

for 

the 

curve 

c 3 , 

p « 

0.631. 


(IV-4) 


Various values of powers P can be tried in the theoretical method in order 
to find out desirable value for better overall correlation with experimental 
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FIGURE IV-19 - VELOCITY PROFILE SIMILARITY IN THE WAKE IN THE VICINITY OF TRAILING 
EDGE OF SINGLE-COMPONENT AIRFOIL WITH SHARP TRAILING EDGE 









data. The reason for the larger scatter in experimental data for similarity 
curves of Figures I V— 1 9 and IV-20 is that these measurements were performed by 
the use of, total-static pressure probe combination in the region where the circu- 
latory flow may be present. This is due to the mixing of upper and lower surface 
boundary layer and the airfoil trailing edge thickness —even though it is small 
it has a finite value. The measurements by the total-static pressure probe, such 
as shown in Figure I I 1-9, in circulatory flow are not very accurate. 


IV. 2. 2 Generalized Parameters for the Pressure Distribution fn the Wake of 
Single Component Airfoil wTth Sharp TraiTing Edge : F i gu res I V- 2 1 through IV-24 

show the plots of non-dimensional pressure or velocity distribution at upper and 
lower edges of the wake for angles of attack ot of 0.0°, 8.0°, 8.0° (fixed 
transition) and 10.79°, respectively. Parameter for the non-dimensional pressure 
at the edges of the wake is defined as 


Y 



(IV-5) 


where U e = velocity at edges of wake 
U = freestream velocity 

and subscripts T. E. = values at or in the neighborhood of trailing edge 

over airfoil surface, and 
(x) = chordwise locations in the wake. 

Non-dimensional parameter y was chosen because, as seen from equation ( I V- 5 ) , 
its value along the wake will vary from zero at the trailing edge to the value 
of 1.0 at very far distance from the trailing edge. Experimental data, shown in 
Figures IV— 21 , IV-22, 1V-23 and IV-24, indicate that when values of y u or yl are 
plotted against distance X'/C, measured from the a i rfoi 1 t ra i 1 i ng edge, in the 
wake then experimental data for upper wake and lower wake arrange nicely to form 
smooth curves for all angles of attack shown. However, different curves are 
obtained for upper and lower edges of wake and moreover curves of y vs. X'/C are 
different for various angles of attack. Thus, in this form y vs. X'/C relation- 
ship cannot be utilized as an ingredient to the theoretical drag computation 
method. By the use of physical reasoning and dimensional analysis, a trans- 
formed X'/C coordinate for the wake flow was arrived at, such that when y u or 
Y|_ is plotted versus transformed coordinate along the wake, then experimental 
data, for both upper and lower edge of wake and for all angles of attack, 
arranged on a single curve. This dimensionless transformed X-coordinate is 
given by the following equation: 
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FIGURE IV-20 - VELOCITY PROFILE SIMILARITY IN THE WAKE IN THE NEIGHBORHOOD OF 

TRAILING EDGE OF SINGLE-COMPONENT AIRFOIL WITH SHARP TRAILING EDGE 
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where X= transformed X-coordinate along the wake 


X' = distance along wake measured from the airfoil trailing edge 
c = airfoil chord 1 ^ 


<S U = boundary layer thickness at the trailing edge on upper surface 
of ai rfo{ 1 

5|_ = boundary layer thickness at the trailing edge on lower surface 
of airfoil 


6q = greater value of 6 U or 5[_ 

Figure 1 V— 25 shows the plot of dimensionless pressure distribution at the 
edges of the wake versus dimensionless transformed coordinate X. This figure 
shows that all points at different chordwise locations and various angles of 
attack arrange nicely on a single curve. The functional relationships between 
parameters y and X is given by 

= i 0.132 

Yu *” (0.363 + X u ) 2 


Figure 1 V— 26 shows the generalized parametric representation and functional 
relationship for the pressure distribution on the locus of minimum velocity in 
the airfoil wake for sharp trailing edge airfoil. The parameters for this uni- 
versal pressure distribution along the locus of minimum velocity point were 
arrived from the consideration of the flow behind the backward facing step and 
from physical reasonings. In this case also, by the choice of proper trans- 
formed X-coordinate, experimental points for Cp,j m . n arrange themselves very, 
nicely on a single curve. The functional relationship between the parameters 
for the pressure along the locus of minimum velocity in the airfoil wake behind 
the sharp trailing edge airfoil is given by. 


- V Um - j - n -(*) » 0.556 e"* 162$ + (0.444 + 1 .121 B + 0.2853 2 )e~ 23 (IV-8) 

C P T.E. * 


where 


Cpu min (x) 


= static pressure coefficient along the locus of minimum 
velocity in wake 


Cp j g = pressure coefficient at the trailing edge of airfoil 

£ = transformed X coordinate for Cp along locus of minimum 
veloci ty 
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FIGURE IV-25 - DIMENSIONLESS UNIVERSAL P 
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g* 1 T - = sum of displacement thickness on upper and lower surface 

° a ’ at the trailing edge of the airfoil 

x = distance along chordline in wake from airfoil trailing edge 

c = a I rfoi 1 chord. 


IV. 2. 3 Generalized Parametric Representat ion of Shear Stress at Velocity 
Minimum :. The theoretical Integral equations of Section If for the solution of 
viscous flow in airfoil wake, contained coupling terms involving shear stress at 
the velocity minimum. It was noted in previous paragraphs that velocity profiles 
in the airfoil wake were ‘‘similar 11 for proper choice of similarity function and 
similarity variable. This fact suggests that it is possible to arrive at 
generalized parameter for the representat ion of shear at minimum velocity point 
from Prandtl's new shear stress hypothesis for free jet flows. According to this 
hypothesis, at any X-location the turbulent shear stress is given by, 

T (y) = K lPc b C U C ^ ( I V— 9) 

where x(y) = shear stress at any y ordinate in viscous layer 
b c = characteristic width of viscous layer 
U c = characteristic velocity in viscous layer 
K x = constant. 

For the upper wake the suitable characteristic width is (nomenclature is 
shown in Figure 1V-27), 


b c - (Y 4 - Y 3 ) (x) 

and the characteristic velocity is 

U ° = Ueu (x) ~ Umin (x) 

For the lower wake, the characteristic width is 

b - (Y S - W M 

and the characteristic velocity is 


u c = Ue L 


U„ 


(x) mm (x) 


(IV-10) 


(IV-11) 


(IV-12) 


(IV-13) 
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The locus of the minimum velocity U m y n , , is contained in both upper wake 
and lower wake layers. Assuming that shear^ x at the locus of minimum velocity is 
proportional to the product of the four quantities given by equations (lV-10), 

( IV— 11)* (IV-12) and (IV-13)» the following functional representation can be 
written after non-dimensional i zing, 


T(y 3 ) 

ipU e 2 


~ Y3) (Y3~Y2) i U e U ( x ) U min( x ) Ue L( x ) U min( x ) 

<Y4 - y 2 ) * (Y4 - y 2 ) ** u, r\ u ro " 


(IV-14) 


Figure IV-27 shows the functional relationship between the non-dimensional 
shear on the locus of the minimum velocity and the product of the parameter indi- 
cated in equation (IV-14). Experimental data plotted in this figure are obtained 
during the present studies in the wake of sharp trailing edge test airfoil at 
various chordwise locations and various angles of attack. Non-dimensional shear 
t ( y 3 )/(ipU e 2 ) was obtained from the indirect shear measurements from measured 
experimental velocity profiles and pressure distributions by the aid of the 
numerical method. The curve fit for the functional relationship of Figure IV-27 
is given by 


T (y 3 ) 

iPUe 2 


115.2 e 2 


(iv-15) 


where T(y ) = shear stress on the locus of minimum velocity 


and 


= arithmetic mean of velocities at upper and lower edges of wake 


r (Y4 ~ Ya) (Y3-y 2 > | U e u ( x ) u min( x )w U e L ( x ) U ™in( x ) ) ~l 

|_(y4-y 2 ) (y 4 - y 2 ) 1 n fj 


IV . 3 Presentation of Correlation Results 

Computer program subroutines were developed for the purpose of comparing 
the results of theoretical computation method, developed under this study, with 
the experimental data for three airfoil configurations. Input to the computer 
program consists of initial conditions of boundary layer and pressure quantities 
on upper and lower surfaces of the airfoil in the vicinity of the trailing edge. 
Experimentally measured values for these input conditions were used as input to 
the computer program. Computations for the characteristics of the flow in the 
wake and profile drag were performed at several angles of attack for the 
following three airfoils: (i) Joukowski airfoil with thickness ratio t/c of 12 

percent, (ii) present sharp trailing edge test airfoil with thickness ratio of 
15 percent, and (iii) NACA 63 1 -01 2 airfoil. Results of correlations are 
discussed in the following paragraphs. 
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(i) Joukowsk i Ai rfo i 1 : Figure IV-28 shows the results of computations for 

upper and lower edge of velocity profiles in the wake of Joukowski airfoil at 
zero degree angle of attack. In this figure the locus of minimum velocity, 
which in this case is a straight line because of symmetrical a i rfoi 1 at zero 
angle of attack, is also shown. Experimental data shown plotted in this figure 
by symbols indicate that good correlation is obtained. 

; • 4 - 

Figure IV-29< shows the plot of experimental measurements (Reference 7) of 
boundary layer momentum thickness, displacement thickness and the form factor 
a£ the trailing edge of Joukowski 's a i rfoil'i a s v ^funct ion 6 fan angle of attack. 
This plot shows that at an angle of attack of 6°,' the lower surface boundary 
layer is separated near the trailing edge as indicated 6y the value of form 
factor H] ower = 2.15. At an angle of 9°, both upper and lower surface boundary 

layer are separated; this is indicated by the values of form factors and also 

by sudden increase in slope d<5*/da for ^the upper surface at a = 9°* These 
values of physical boundary layer thickness and pressure coefficient at the 
trailing edge were input to the profile drag computer program. Values of 
boundary layer physical thicknesses were computed from experimental momentum 
thicknesses and form factors by the use of the following equation: 

s 1 + 0.1 M e 2 

? = 0.17 Ti - exp {- 3 » 5 ( H- 1 ) }] (IV-16 

Figure IV“30 shows the computed variations of displacement thickness, 
momentum thickness and form factor in the wake of Joukowski airfoil at an angle 

of attack a = 0°. Figure IV-31 shows the plot of above parameters at an angle 

of attack a = 6°. Experimentally measured values of these parameters, which 
have been shown plotted in Figure 1V-30 and IV-31 as symbols, indicate that 
good agreement is obtained between computational results and experimental data. 
Figures IV -32 and I V- 33 show the comparisons between results of comptations of 
velocity profiles and experimental data in the wake of Joukowski airfoil for 
angles of attack of 0° and 6°, respectively. The above comparison is made at 
X’/C = 0.1, 0.25 and 0.5 in Joukowski airfoil wake and this comparison between 
computed and experimental velocity profiles indicate that reasonable correlation 
is obtained. ' 

Figure IV-3*f shows the plot of computed values of the profile drag coef- 
ficients as a function of an angle of attack and comparison with experimental 
data. This figure shows that the comparison between computed profile drag 
values and experimental data is reasonable at angles of attack a = 0°, 3° and 
6° whereas at a = 9° the comparison between present theoretical computational 
results and experimental data is not so good. Experimental measurements at 
airfoil trailing edge, which is shown in Figure IV-29, indicate that at a - 0° 
and 3° flow separation is absent on airfoil surfaces. Incipient separation 
exists on a i rfoi 1 surfaces near the trailing edge in the neighborhood of a = 6° 
and appreciable flow separation probably exists on airfoil surfaces for angles 
of attack greater than a - 6° as evidenced from Figure IV-29* The discrepancy 
in "drag between results of computations and experimental drag data for angles 
of attack greater than 6° can be due to the following two reasons. The initial 
conditions of boundary layer quantities, which are required as input to 
computer program were obtained from data of Figure IV-29- These data, which 
were obtained from Reference 7, were measured by pressure probe and the validity 
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FIGURE IV-34 “ COMPUTED VARIATION OF PROFILE DRAG COEFFICIENTS WITH ANGLE OF ATTACK 
FOR JOUKOVSKI AIRFOIL AND COMPARISON WITH EXPERIMENTAL MEASUREMENTS 



of measurements by the pressure probe In the separated boundary layer is 
questionable. Secondly, the validity of assumptions which were used in deriving 
theoretical wake flow equations of Section II, when appreciable flow separation 
exists on airfoil surfaces, is not known at the present time. For the purpose 
of clarifying the above-mentioned situation, need exists for (i) developing 
proper experimental techniques for obtaining accurate and reliable measurements 
in separated flow region and (ii) then making appropriate modifications to the 
theory for wake flow calculations as are necessary for accurate drag predictions 
in the presence of flow separation. 

(ii) 'S harp T.E. Test Airfoil : Figure IV-35 shows the plot of computed 

characteristic loci at an angle of attack a = 8° and for the case of free tran- 
sition on the surfaces of the test airfoil. Computed loci are plotted for upper 
edge of the wake, lower edge of the wake and the locus of minimum velocity for 
the velocity profiles in the test airfoil wake. Experimental values, which are 
shown in this figure as symbols, indicate that agreement between computed values 
for these characteristic points in velocity profiles and experimental data is 
good. It should be noted that slope of these loci downstream of X'/C of approxi- 
mately 0.5 become approximately parallel to the free stream direction. 

Figures IV- 36 , IV-37 and IV -38 show the plots of computed values of momentum 
thickness and form factor distribution in the wake of the test airfoil at angles 
of attack a = 0°, 8° (fixed transition) and 10.79°> respectively. Computed 
values of momentum thickness distribution is shown plotted for lower half of wake 
profile and total wake velocity profile'; difference between these two momentum 
thickness distributions is the upper wake momentum thickness. Computational 
results for form factor distribution in airfoil wake indicate that even though 
the form factor on the upper and lower surfaces at the trailing edge are vastly 
different, the values of form factors for upper and lower half of wake velocity 
profile become identical at a very short distance from the airfoil trailing edge. 
For this reason only one value of computed form factor distribution for the 
flow in the wake is shown in Figures IV- 36 , IV-37 and IV- 38 . Experimental data 
for the form factor distributions, which are shown plotted in this figure, are 
either for lower half wake or upper half wake because measured values of form 
factors are very nearly equal for either upper or lower wake velocity profile. 
Figure I V- 39 shows the plot of results of calculations of profile drag versus 
angle of attack for the test airfoil. Experimental data which are shown plotted 
in this figure indicate that agreement between experiment and theoretical method 
is reasonable. 

(iii) NAC 63 1 “ 01 2 Airfoil - Results of correlation between computations 
and experimental data tor this airfoil are shown in Figure IV-40. Experimental 
data shown in this figure are obtained from Reference 5- For the purpose of 
correlation for this airfoil, measured values of boundary layer quantities on 
airfoil upper surface at trailing edge were used as input to the profile drag 
computer program. Boundary layer quantities on the lower surface at the 
trailing edge of the airfoil which are also required input to the computer 
program, were obtained from the output of the computer program subroutine for 
ordinary boundary layer calculation method. Experimental pressure distributions 
on the lower surface of the airfoils were used as an input for the purpose of 
calculation of lower surface boundary layer quantities at several angles of 
attack. In addition, values of experimental pressure at the trailing edge of 
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FIGURE IV -36 - COMPUTED VARIATION OF 0/C AND H IN THE WAKE OF SHARP T.E. TEST 
AIRFOIL @ K - 0°AND COMPARISON WITH EXPERIMENTAL MEASUREMENTS 








FIGURE IV-37 - COMPUTED DISTRIBUTION OF ©/C AND H IN THE WAKE OF SHARP T.E. TEST AIRFOIL 
AT e*= 8.0° (FIXED TRANSITION) AND COMPARISON WITH EXPERIMENTAL DATA 
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the airfoil were used as input to the profile drag computer program subroutines. 
Figure I V— 40 shows the comparison of the results of computations for the profile 
drag coefficient of NACA 63 i -01 2 airfoil with experimental data. Plot of 
experimental C|_ vs. a for this airfoil is also shown in this figure; in addition, 
values of form factor on the upper surface of the airfoil at its trailing edge 
are noted on this Cl - a curve at several angles of attack. This is done for the 
purpose of indication of approaching trailing edge separation and also for the 
purpose of indicating values of angles of attack for which separated flow condi- 
tions exist at airfoil trailing edge. The comparison of computed profile drag 
quantities with experimental data is good as seen from this figure. 


V. CONCLUSIONS AND RECOMMENDATIONS 


From the theoretical and experimental studies presented in this report it 
is possible to make the following conclusions and recommendations for the future 
studies. 


V. 1 Concl us ions 


(1) The total-static pressure probe combination, such as used in the 
present studies, can be used for reliable measurements in viscous flow over air- 
foil surface when this flow is attached or non-separated. However, the measure- 
ments by such a probe, for the flow approaching separation or in separated flow 
region, are subject to interpretation and hence can be used only in qualitiative 
sense. 

(2) Measurements of velocity profile by the use of hot-wire anemometer 
technique have been found more reliable, both qualitatively and quant i at i vel y , 
than similar measurements performed by the use of pressure probe. However, with 
the present state of the art in hot-wire anemometry technique, measurements of 
velocity profiles in the region of low values of velocities (less than 30 fps 
and nominal values of turbulence level) are not very accurate. As the hot-wire 
anemometer measures only the absolute magnitude of velocity then with the above- 
mentioned present limitation, it is difficult to determine the separated flow 
velocity profile with the desired accuracy. 

(3) Shear stress profiles for the wake flow have been determined by the 
use of indirect measurements and these values are used in the parametric form 
in the present theoretical method. The accuracy of indirect shear measurements 
for the wake flow in the trailing edge region is questionable because of the 
existence of the circulatory shear flow in this region. For this reason it is 
recommended to perform the direct shear measurements by the use of such devices 
as hot-wire X probe or by the use of Laser Doppler Velocimeter and make 
comparisons with the present indirect shear measurements. 
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(A) From the comparison between results of computer program subroutines 
with experimental data on three airfoil configurations, the following three 
remarks can be made: 

(i) Agreement between computed and experimental velocity profiles for 
the flow in the wake is quite good for symmetrical wake at large distances from 
the trailing edge. However, for the unsymmetrical wake and especially in the 
region of wake near the airfoil trailing edge, the discrepancies exist between 
the results of computations and experimental measurements. These discrepancies 
might be due to inaccuracies in measurements or theoretical calculations or 
both in this region. 

(ii) The computed variations of integral quantities in the airfoil 
wake, such as momentum thickness, displacement thickness and the form factor, 
agree well with experimental data for distances far from the trailing edge. In 
the region near the trailing edge, however, discrepancies are observed between 
the theory and experiments. 

(iii) Results of computations of profile drag agree quite well with 
experimental measurements in the range of values of angles of attack from 

a = 0° to a corresponding to the occurrence of incipient separation on the 
upper surface of the airfoil. However, for high values of angles of attack 
where large region of flow separation is present on the airfoil surface, com- 
puted values of the profile drag are not in good agreement with experimental 
measurements. This discrepancy can be attributed to the breakdown of the 
assumptions, in case of separated flow, which were used in deriving theoretical 
equations for the present method. 


V. 2 Recommendations 


(1) For the purpose of establishing the limitations and restrictions of 
the present method, it is necessary to perform correlations with experimental 
data for the wake flow and profile drag on several classes of airfoil configu- 
rations. Thus performing calculations on supercritical and conventional 
airfoils of several thickness ratios and camber distributions and at several 
values of lift coefficients and comparison of computational results with 
experimental measurements would be very valuable for the purposes of needed 
refinements for this new method. 

(2) This method is developed for sharp trailing edge airfoils and for 
low subsonic free stream Mach number where the flow can be considered essen- 
tially incompressible. However, for the general applicability in the design and 
development work on wing sections it is recommended that the present method be 
extended for high subsonic flows and for airfoils with thick trailing edges. 

(3) Experimental and theoretical studies are recommended for the refine- 
ments of the present method for the purpose of the prediction of profile drag 
when extensive region of flow separation is present on airfoil surface. Such 
conditions exist at Cl^x conc *itions on airfoils exhibiting trailing edge stall. 
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(4) The basic approach of the theory of the present method is valid for 
the more important situations of the computation of profile drag for multi- 
component airfoil sections. Under the present study the validity of this 
approach is proven for single-component airfoil sections. It is hence recom- 
mended that this approach be extended to airfoils with more than one component. 
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APPENDIX A 


DETERMINATION OF SIDE WALL BLOWING REQUIREMENTS 


It is known that in the case of two-dimensional testing at conditions of 
high lift, the validity and accuracy of experimental measurements is sometimes 
questionable because of the existence of boundary layer separation on the test 
section sidewalls. The presence of the two-dimensional model at moderate and 
high angles of attack in the wind tunnel test section creates high adverse 
pressure gradients in the direction of flow on the side walls of the tunnel test 
section upstream of the leading edge of the model. This high adverse pressure 
gradient causes the separation of boundary layer on the wind tunnel test section 
side walls at the juncture of the model with the side walls. Separated boundary 
layer on the side walls then flows toward the center of the model with the 
result that the flow pattern on the airfoil model is three-dimensional thus 
making it impossible to obtain measurements of two-dimensional aerodynamic 
characteristics of given airfoil model at moderate and high angles of attack. 

The above-mentioned phenomena has been observed by various investigators at NASA, 
R.A.E. in England, and at Lockheed- Georg i a Company. 

The solution of the above problem can be approached by the use of several 
methods, such as distributed blowing through slits located in the sidewalls at 
juncture of the model with side walls, distributed suction through pores or 
slits located upstream of the model or at the juncture of the model with the 
tunnel side walls, and blowing high energy air through a single slit located at 
an appropriate distance upstream of the leading edge of the model on both side 
walls. For the present investigations boundary layer control on wind tunnel 
side walls was accomplished by blowing the appropriate amount of high energy air 
through single slits located at approximately one chord length upstream of model 
leading edge on both side walls. Figure MI-8 shows the general layout of the 
side wall boundary layer control system used in the present investigation. The 
following paragraphs describe the procedure that was used in determining the 
appropriate amount of high energy air for various model configurations and 
angles of attack. 

Velocity profiles on the side walls of the tunnel were measured without the 
presence of the model at locations corresponding to the approximate locations of 
the trailing edge of the main component and flap. Then, with the model in the 
tunnel in various configurations and angles of attack, measurements of velocity 
profiles were made on the side walls of the tunnel at the same locations. These 
measurements were used to determine blowing requirements. It is necessary that 
the side wall BLC be just sufficient to prevent boundary layer separation on the 
side walls. It has to be emphasized that the amount of blowing on the wind 
tunnel side walls for boundary layer control has to be precise — for example, if 
there is excess blowing then the value of freestream q^ in the test section 
would be affected. In addition, an excess amount of blowing on wind tunnel walls 
creates an undesirable pressure field around the two-dimensional model; this is 
due to the fact that vortices are created due to the rubbing action of low 
velocity tunnel freestream air with the high velocity blown air. These vortices 
give rise to tip effects or the effects of finite span on two-dimensional test 
model. If the amount of blown high energy air on the tunnel side walls is 



insufficient then, of course, the boundary layer on the side walls separates. 

In order to avoid both undesirable effects, the proper value of the pressure 
ratio P ts , P Ss and the proper geometrical dimensions of the slits are required. 

A simplified theoretical equation of the amount of blown air required to 
suppress boundary layer separation on wind tunnel side walls can be derived as 
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where A0 


max 
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h s 

Ps 
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maximum value of the difference in momentum thickness on wind 
tunnel side walls in the presence of and absence of model in the 
test section, inches 

airfoil chord, inches 

effectivity factor 

height of blowing slits, inches 

density of high energy air in blowing plenum 

freestream density. 


The blowing momentum coefficient C in equation (A-l) is defined as 




(mass of blown air) * (velocity of blown air) 
1 2 

y p U (wing chord) • (span of model) 


(A-2) 


The value for the ratio of velocity of blown air through the blowing slit 
to the freestream velocity, V s /U , is approximately by 


D 0.286 -i 1/2 




where Y = ratio of specific heats for air (1.4) 

R = gas constant for air (=53*3) 
g c = gravitational constant (= 32 . 2 ) 

Ps s = static pressure at the exit of blowing slit, psia 
PT S = plenum air total pressure. 


(A-3) 


The amoung of that is available corresponding to a given pressure ratio, 
Pt s /Ps s , is given by the following equation: 
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(A-4) 


where Vg/U^ is given by equation (A— 3 ) 
ht = span of the model , inches 

= 30 inches for the present tunnel 
Kw = width of blowing slit, inches 
= 30 inches for the present tunnel 

Figure A-1 shows the plot of available C y for the unit slot height as a 
function of pressure ratio Pt^/PSg. Figure A-2 shows the plot of velocity ratio 
as a function of pressure ratio Pt s /Ps s . Figure A~3 shows the plot of effectivity 
factor as a function of velocity ratio U ei /V s . The curve of Figure A-3 is ob- 
tained from knowledge of experimental measurements and is empirical. The curves 
shown in Figures A-1 and A-2 are plots of equations (A-4) and (A-3), respectively; 
these curves are constructed for freest ream dynamic head q^ of 60 psf and free- 
stream total temperature of 500°R. By simultaneous use of equation (A-1) and the 
curves of Figures A-1, A-2, and A~3, a desired value of the pressure ratio of 
blown high energy air can be calculated by trial and error such that the required 
Cy is equal to the Cy which is available corresponding to a given pressure ratio 

Pt s /Ps s • 
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APPENDIX B 


ADDITIONAL INSTRUMENTATION 


The test section flow conditions and wind tunnel performance data were 
monitored with conventional instrumentation. Fan speed, clutch current, and 
tunnel airstream total temperature were displayed on the control console near 
the speed control potentiometer as shown in Figure B-l . During these tests, 
atmospheric pressure was measured using a 9-inch diameter Wallace and Tieman 
absolute pressure gage with reading increments to 0.01 pounds per square inch. 
Atmospheric pressure was manually input to the data system before each run. 

The test section dynamic pressure was set and monitored visually with a 70- 
inch water manometer connected to pressure orifices located upstream and down- 
stream of the wind tunnel contraction. The manometer reading gave the static 
pressure drop across the contraction. This reading was related by calibration 
to the dynamic pressure in the empty test section. The manometer, which is 
shown in Figure B-2 has provisions for reading the height of the water column 
to within ±0.001 inch with good repeatability. This manometer was also used to 
calibrate some of the pressure transducers used during these tests. 

The contraction pressures were also read in electrical units using two 
Statham type PM6TC pressure transducers rated at ±1.0 pound per square inch 
differential. These transducers can be expected to be accurate to within about 
0.5 pounds per square foot in the normal operating range. One transducer was 
connected to the contraction pressure orifices in parallel with the water manom- 
eter and was calibrated versus dynamic pressure in the empty test section. The 
other was connected to the upstream or high pressure end of the contraction and 
was calibrated versus test section total pressure. The output of these trans- 
ducers was transmitted to the data acquisition system for recording and use in 
data reduction. 

The wind tunnel airstream total temperature is measured by a thermistor 
mounted on a probe in the settling chamber. Temperature information from the 
probe is displayed on the control console and is transmitted electrically to the 

data acquisition system. A diagram of the basic wind tunnel instrumentation is 

shown in Figure B-3. 

Static pressures were measured at 41 orifice locations on the model using 
two Statham model D3~GM scanivalves with two PM -131 pressure transducers. The 
scanivalves were ganged together and actuated by a single solenoid type stepper. 

Valve number one contained a transducer rated at 12.5 pounds per square inch 

and was connected to the orifices on the model upper surface. Valve number two 
contained a 2.5 pounds per square inch transducer and was connected to the 
orifices on the lower surface of the model. The output from these transducers 
was recorded by the data acquisition system for use in computing pressure co- 
efficients, lift coefficients, and local velocities. Figure B-*t shows the 
scani valve instal 1 at ion. 

Total pressure in the boundary layer control plenums was calibrated versus 
the pressure in the auxiliary air supply duct between the plenums and the control 



valve. The supply duct pressure was displayed on an ordinary 0 to 100 pounds 
per square inch pressure gage mounted on the control console as shown in Figure 

B-1. 
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APPENDIX C 

DATA ACQUISITION SYSTEM 


General Description 


*The data processing system utilized for this contract was set up for two 
specific purposes, these being the measurement , of the static pressure distribu- 
tion on the surface of the model and the measurement of the velocity profiles 
within the boundary layer and wake of the model. The measurements were made 
using a Data Acquisition Unit (D.A.U.) which was controlled by a real-time 
digital computer which activated scanivalve units for obtaining the static 
pressure distribution and traversed a pressure or hot wire anemometer probe for 
the velocity profiles. 

The system used for acquiring and reducing the test data is shown in 
Figure C-1. The heart of the system is a Lockheed Electronics MAC 16 computer. 
The raw data was made available, in an abbreviated form, on a teletype for on- 
line monitoring of the test, and in its entirety on paper tape using a high 
speed punch. Final data reduction was accomplished on a UNIVAC 1106 Central 
Computing System with remote access terminals. 

An additional MAC 16 computer was used in an executive mode to permit time 
sharing of the main computer. The disc storage associated with this computer 
was intended to provide the basis of an on-line data reduction capability in 
conjunction with a UNIVAC 4 1 8 . During this contract, however, the immediate 
availability of the remote access terminals made the on-line capability 
unnecessary. 


The Data Acquisition Unit 

The data acquisition unit incorporated an analog multiplexer 1 for transmis- 
sion of measured data and a digital multiplexer for transmission of thumb-wheel 
switch inputs. One additional analog channel was provided for probe position 
measurement. 

A rotary switch permitted any of the measured data to be monitored on a 
digital voltmeter during the test. Variable gain amplifiers in all of the 
analog channels enabled gain levels to be selected that approximated the raw 
data to engineering units for ease of monitoring of the digital voltmeter and 
on-line teletype output. 

The tunnel conditions of 1 q 1 , 'H 1 and temperature were acquired at each 
scanivalve position during the static pressure distribution and at each probe 
position during the boundary layer survey. The thumb-wheel switch inputs 
included a probe option for choice of no probe, pressure probe or anemometer 
probe and also a selection of the number of scanivalves and ports to be cycled 
through the static pressure distribution. Any selection below the maximum 
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required (2 scanivalves and 36 ports) inhibited printout of the pressure data 
during reduction with significant saving in reduction time. This option was 
exercised when a series of velocity surveys were made with no change of basic 
configuration, and hence pressure survey, in which case the minimum number of 
ports was selected that would include the static pressure port at the chordwise 
location of the velocity survey. 

Control of the data acquisition was achieved with one analog channel which 
activated the probe positioning servo mechanism and one digital channel which 
stepped and homed the scanivalves and controlled the multiplexing. 
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